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V 

Heat  transfer  effects  of  longitudinal  vortices  embedded  within  film: 
cooled  turbulent  boundary  layers  on  a  flat  plate  were  examined  for 
freestream  velocities  of  10  m/s  and  15  m/s,  and  for  blowing  ratios  ranging 
from  0.47  to  1.26.  Moderate  strength  vortices  were  employed  having  circu¬ 
lation  to  freestream  velocity  ratios  of  about  1.6  cm.  Spatially  resolved  heat 
transfer  measurements  from  a  constant  heat  flux  surface  and  mean  tempe¬ 
rature  distributions  in  snajiwise  planes  show  that  local  heat  transfer  is 
significantly  affected  by  spanwise  vortex  position,  and  blowing  ratio.  Of 
particular  significance  are  boundary  layer  and  vortex  structural  changes 
which  occur  at  high  blowing  ratios,  fo  ' '  1 
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I.  INTRODUCTION 


The  increasing  need  for  greater  efficiency  in  gas  turbine  engines  has  re¬ 
sulted  in  higher  turbine  inlet  temperatures.  Consequently,  combustors  and 
turbine  blading  are  subject  to  greater  amounts  of  thermal  stress,  thermal  fa¬ 
tigue.  and  creep.  At  present,  gas  turbines  such  as  those  associated  with  mi¬ 
litary  applications,  have  inlet  temperatures  as  high  as  1800-2000  ^C,  [Kef  1 1 
While  the  development  of  improved  alloys  with  higher  melting  points  is  part 
of  the  solution,  the  development  of  efficient  cooling  systems  is  just  as 
important  [Ref.  2j.  In  order  to  design  an  efficient  cooling  configuration,  the 
heat  transfer  distributions  for  the  gas  turbine  components  are  needed 
Because  of  the  complex  geometries  and  flows  involved  near  blades  and  end- 
walls,  accurate  convective  heat  transfer  rates  are  difficult  to  obtain.  [Ref.  31 

Ongoren,  [Ref.  4],  described  the  flow  in  a  turbine  cascade.  As  the  inlet 
boundary  layer  approaches  the  blade,  just  in  front  of  the  blade,  a  horseshoe 
vortex  forms.  At  the  saddle  point,  it  splits  into  a  vortex  on  the  suction  side, 
and  a  vortex  on  the  pressure  side.  The  pressure  side  vortex  becomes  the 
passage  vortex,  moving  from  the  leading  edge  of  the  blade  towards  the  low 
pressure  side  of  the  adjacent  blade.  As  the  suction  side  vortex  convects 
along  the  blade,  it  is  eventually  pushed  away  by  the  passage  vortex  from  the 
adjacent  blade.  A  smaller,  corner  vortex  rotates  in  opposite  direction  to  the 
passage  vortex  as  was  verified  by  Sieverding,  [Ref.  3). 

The  passage  vortex  is  an  example  of  a  longitudinal  vortex  em  edded  in 
a  lilm  cooled  turbulent  boundary  layer.  In  most  cases,  an  embedded  vortex 


WA'  V  *w"  V'  *-  v-"  V<  'O  x-'  \WW-  <w  »V.TV  T»  T«  w\  v,  y .-V.-  «V  YVV. 


ai r.VL'Tft.v'vvwLr^r^n^  rk*-j»  ^  ^ 


has  a  cross-stream  length  scale  approximately  equal  to  the  boundary  layer 
thickness.  Therefore,  the  vortex  is  capable  of  strongly  perturbing  the 
boundary  layer  structure  and  modifying  the  heat  transfer  characteristics.  In 
addition,  longitudinal  vortices  usually  maintain  their  coherence  over  a  long 
streamwise  distance,  meaning  that  the  heat  transfer  effects  behind  an  effec¬ 
tive  vortex  generator  are  likely  to  be  very  persistent.  (Ref.  61 

Film  cooling  is  used  to  provide  a  layer  of  cool  fluid  between  a  surface 
and  high  temperature  free  stream  gases  to  which  it  is  exposed.  The  film 
cooling  not  only  acts  as  a  thermal  insulator  but  also  as  heat  sink  lor  the  hot 
freestream  gases.  The  overall  effect  of  the  film  cooling  is  to  reduce  the 
temperature  of  the  developing  boundary  layer,  which  in  turn  reduces  the 
heat  transfer  to  the  surface.  In  the  present  study,  the  heat  flux  is  calculated 
from  : 

q  ■  h(Tv  ■  Tpoo)  (eqn.  1  1 ) 

where  h,  includes  the  effect  of  the  film  cooling  for  constant  wall  heat  flux,  q ", 
and  constant  free  stream  temperature,  Troo.  In  the  present  experiment 
variations  of  h  correspond  to  variations  of  the  wall  temperature,  Tw.  When 
T^  *  Tfoo,  equation  1.1  indicates  that  q"  -  0,  which  may  not  be  the  case  when 
film  cooling  is  present. 

Goldstein  and  Chen,  (Ref.  21,  defined  an  adiabatic  wall  effectiveness  for 


film  cooling  using: 


The  heal  transfer  is  then  calculated  from  : 

q  =  hf(  Tv  -  Tav)  (eqn.  I  3 J 

Here,  q  "=  0  when  Tv  =  Taw  with  film  cooling,  Without  film  cooling  Tcoo  »  Tav, 
and  equations  1.1  and  1.3  are  the  same. 

Numerous  studies  have  been  conducted  on  the  effects  of  film  cooimg  on 
heat  transfer,  effects  of  secondary  flows  on  heat  transfer,  and  more  recently, 
on  effects  of  film  cooling  and  secondary  flows  on  heat  transfer  in  a  turbuient 
boundary  layer.  Summarizing  the  results  of  experimental  research  on  the 
effects  of  secondary  flows  in  turbine  blades  passages  over  the  past  decade, 
Sieverding  (1984)  concluded  that  it  is  absolutely  essential  to  establish  the 
significance  of  secondary  flows  at  off  design  conditions,  in  particular  in 
regard  to  leading  edge  vortices. 

Studying  the  effects  of  film  cooling  and  secondary  flows  on  heat  trans¬ 
fer,  Golstein  and  Chen  (1985)  performed  an  experimental  study  on  the  influ¬ 
ence  of  the  end  wall  on  film  cooling  of  gas  turbine  blades  using  a  single  row 
of  injection  holes.  Two  years  later,  the  same  authors  performed  a  similar 
study  but  this  time  employed  two  rows  of  injection  holes  for  the  film  cooling, 
[Ref  7J.  From  these  two  studies,  the  authors  concluded  that,  in  the  convex 
side  of  the  blade  there  is  a  triangular  region  where  coolant  is  swept  away 
from  the  surface  by  the  passage  vortex,  while  the  concave  side  was  not  sig¬ 
nificantly  affected  by  secondary  flows  originating  near  the  endwall.  Other 
study  on  film  cooling  effects  is  given  in  [Ref.  8). 


Eibeck  and  Eaton  (1987),  conducted  a  study  of  a  longitudinal  vortex  em¬ 
bedded  in  a  turbulent  boundary  layer  over  a  constant  heat  flux  plate.  The 
authors  found  significant  increases  and  decreases  in  local  Stanton  numbers, 
due  to  the  thinning  of  the  boundary  layer  on  the  downwash  side  of  the 
vortex  and  thickening  on  the  upwash  side  of  the  vortex.  S panwise  heat 
transfer  variations  became  larger  as  the  circulation  of  embedded  vortices 
increased.  [Ref.  61 

Investigating  the  interaction  between  a  single  weak  streamwise  vortex 
and  a  two-dimensional  turbulent  boundary  layer,  Russel,  Pauley  and  Eaton 
(1987),  observed  a  rapid  growth  of  the  vortex  core,  and  a  flattening  of  the 
core  shape  when  the  dimension  of  the  core  radius  became  comparable  to  the 
distance  of  the  vortex  center  to  the  surface.  Adverse  pressure  gradient 
caused  an  increase  in  the  rate  of  core  growth,  and  a  stronger  distortion  of  the 
core  shape.  These  authors  also  provide  an  extensive  review  of  work  on 
vortices  in  boundary  layers.  [Ref  91 

In  [Ref.  10],  a  study  was  conducted  on  heat  transfer  distributions  and 
film  cooling  effectiveness  on  the  endwall  of  an  airfoil,  using  a  full  annular 
low  aspect  ratio  vane  cascade.  The  authors  demonstrated  the  importance  of 
the  horseshoe  vortices  and  secondary  flows  on  the  heat  transfer  and  film 
cooling  distributions. 

Recently,  Joseph  (1986),  conducted  a  study  of  the  effects  of  embedded 
vortices  on  heat  transfer  in  film  cooled  turbulent  boundary  layers.  The 
author  showed  that  the  effects  of  the  vortei  on  heat  transfer  are  significant 
and  important :  on  the  downwash  side  of  the  vortei,  heat  transfer  is  aug¬ 
mented,  effects  of  the  film  cooling  are  negated  and  local  hot  spots  will  exist. 
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Near  the  upwash  side  of  the  vortei,  coolant  is  pushed  to  the  side,  appearing 
to  augment  the  protection  provided  by  the  film  cooling,  (Ref.  31.  Evans 
(1987),  studied  the  fluid  mechanics  of  vortices  in  boundary  layers  with  film 
cooling,  and  showed  that  the  vortices  completely  dominate  the  flow  field, 
especially  near  the  downwash  side  (Ref.  1 ). 

The  objective  of  this  thesis  is  to  increase  the  understanding  of  the  effects 
on  heat  transfer,  of  a  longitudinal  vortex  embedded  in  a  film  cooled  turbu¬ 
lent  boundary  layer  as  affected  by  :  1 )  spanwise  vortex  position  relative  to 
the  location  of  film  cooling  injection  holes  and,  2)  blowing  ratio.  These 
effects  are  important  regarding  turbine  blade  and  endwall  heat  transfer. 
The  study  was  carried  out  under  a  series  of  steps.  The  first,  was  the  design 
and  construction  of  the  heat  transfer  test  surface,  followed  by  qualification 
test  to  verify  uniform  wall  heat  flux,  and  an  energy  balance  to  identify  and 
quantify  the  heat  losses.  After  completion  of  the  experimental  apparatus, 
four  types  of  heat  transfer  test  were  conducted  :  ( 1 )  heat  transfer  data  with 
developing  boundary  layer  only,  (2)  boundary  layer  and  embedded  vortei, 
(3)  boundary  layer  with  film  cooling  only,  and  (4)  boundary  layer  with  film 
cocling  and  embedded  vortex. 


II.  EXPERIMENTAL  APPARATUS 

A.  WIND  TUNNEL 

The  wind  tunnel  pictured  in  Figure  3  ,  built  by  Aerolab,  was  designed  to 
provide  a  flow  field  from  the  nozzle  with  uniform  velocity  and  low  turbu¬ 
lence  intensitv.  It  is  designated  the  NPS  Shear  Layer  Research  FjcimI'. 
tSLRFi. 

I .  Description. 

The  SLRF  is  a  wind  tunnel  of  the  open  circuit  blower  type  with  lan 
upstream  and  air  entering  the  blower  inlet  from  the  surrounding  room  The 
air  speed  through  the  test  section  can  be  adjusted  from  5  to  40  m/s.  The 
tunnel  frame  has  leveling  screws  to  adjust  the  center  line  of  the  tunnel  to  a 
horizontal  position  The  discharge  of  the  fan  slips  into  the  inlet  of  the  dif¬ 
fuser  with  a  1.6  mm  clearance  to  isolate  vibrations  from  the  fan  to  the  wind 
tunnel  body  The  diffuser  section  contains  a  filter  pack  and  a  nozzle  leading 
to  test  section.  The  test  section  is  a  rectangular  duct,  3  048  m  long  and 
0.6090  m  wide  It  is  designed  with  numerous  pressure  tabs  and  four  38  x 
20.3  cm  access  ports  along  each  of  the  side  walls  to  provide  easy  access  The 
top  wall  is  a  continuous  panel  fabricated  from  4  76  mm  thick  Lexan  sheet, 
continuously  sealed  with  neoprene  along  the  edges.  The  ceiling  height  is 
adjustable  to  permit  changes  in  the  pressure  gradient  along  the  length  of  the 
test  section.  Additionally,  the  top  wall  contains  numerous  instrument  ports 
for  the  measurements  of  various  flow  characteristics  The  floor  of  the  test 
section  consists  ot  three  0.6096  m  and  one  1  2  1 92  m  long  sections  which  are 


removable  and  replaceable.  These  sections  are  all  0.6090  m  wide  and  are 
sealed  w’lth  O  '  rings  around  the  sides.  Further  discussions  of  the  wrina  tun¬ 
nel  are  contained  in  (Ref.  1 1 1  and  (Ref.  3:  pp.  38], 


A  schematic  of  the  test  section  components  with  the  coordinate 
system  employed  in  the  present  study  are  shown  in  Figure  1.  An  unhealed 
starting  length  oi  1.10  tn  exists  upstream  of  the  heated  test  surface  Che  in¬ 
jection  nozzles  are  located  1.08  m  downstream  of  the  boundary  layer  mo 
and  0.02  m  upstream  of  the  test  surface.  The  leading  edges  of  the  vortex 
generators  are  placed  0.479  m  downstream  of  the  boundary  layer  trip  as 
shown  in  Figure  2. 

2.  Qualification  and  Performance. 

Prior  to  its  relocation,  extensive  qualification  tests  of  the  Shear 
Laver  Research  Facility  were  conducted  by  Ligram,  [Ref.  12].  Results  show 
that  the  variation  of  total  pressure  at  the  exit  plane  of  the  nozzle  is  less  than 
0.4%  at  26  m/s  and  34  m/s.  Mean  velocity  varies  less  than  0.7%  for  the 
same  mean  freestream  speeds.  From  five-hole  pressure  probe  measure¬ 
ments,  the  velocity  angle  deviation  is  nowhere  greater  than  0.6  degrees  at 
the  nozzle  exit  plane. 

Profile  measurements  of  the  mean  velocity  and  longitudinal  turbu¬ 
lence  intensity  in  the  turbulent  boundary  layer  developing  at  20  m/s  indi¬ 
cate  normal,  spanwise  uniform  behavior.  For  this  qualification  test,  and  all 
results  which  follow,  the  boundary  layer  was  tripped  near  the  exit  of  the 
nozzle  with  a  1.5  mm  high  strip  of  tape.  For  the  present  study,  total 
pressure  measurements  along  the  test  section  were  made  after  relocation  of 
the  SLRF,  in  order  to  adjust  the  top  wall  for  zero  pressure  gradient.  The 


maximum  pressure  difference  along  the  test  section  was  0.007  in  of  H?') 
Free  stream  turbulence  intensity  was  measured  to  be  0  00085  1  8  one 

hundredths  of  one  percent  or  0.085 %)  at  20  m/s,  increasing  to  0.00095  at  30 
m/s. 

B  INJECTION  SYSTEM. 

The  injection  system  provides  film  coolant  at  temperatures  above  am¬ 
bient.  The  coolant  is  ejected  from  a  single  row  of  injection  holes  into  the 
boundary  layer  developing  along  the  bottom  wall  of  the  test  section.  The 
diameters  of  the  injection  holes  were  scaled  relative  to  boundary  layer 
thickness  to  be  similar  to  a  turbine  blade,  with  5/d  ratio  ranging  from  0.37  to 
0  40.  The  free  stream  air  is  at  ambient  temperature,  thus  the  direction  of 
heat  transfer  is  opposite  that  of  a  gas  turbine.  The  temperature  difference 
Tw  -  Troo  range  for  this  study  has  been  kept  less  than  20  QC.  to  minimize  ihe 
effects  of  variable  fluid  properties. 

The  injection  parameters  M  and  8  were  scaled  to  resemble  parameters 
near  gas  turbine  blades  where  M  ranges  from  0.47  to  1.26  and  8  varied 
from  1.39  to  1.69.  Since  the  average  temperatures  of  the  coolant,  plate  and 
freestream  were  approximately  the  same  for  different  runs,  Trc  =  51  QC,  and 
Tv  =  40  QC  and  Troo  ~  20  QC,  variation  of  0  was  not  a  parameter  considered  in 
this  study.  Due  to  the  reverse  direction  of  heat  transfer  for  the  present  ex¬ 
perimental  apparatus  Trc :  Troo  :  Tw  ratio  is  1.27  :  0.5- 0.55  :1  0  as  compared 
to  0.67-  0  83  :  1.5  :  10  for  actual  gas  turbines 
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1.  Description 

Air  lor  the  injection  system  originates  in  a  two  stage.  ISO  osig 
tngersol-Rand  air  compressor,  10  HP,  model  *  71TD.  From  the  compressor 
air  Hows  through  an  adjustable  pressure  regulator,  a  cut-off  valve,  a 
reinforce  flexible  tubing  (2.54  cm  inside  diameterf  a  moisture  separator  a 
i low  regulator,  a  fisher  and  Portor  rotometer  (lull  scale  9  3*4 5K--3  m^/s  iv  t 
5oFM.  model  -  1 0.43565  At  a  dil fuser,  and  finally  into  the  injection  iieat 
exchanger  and  plenum  chamber.  The  rotometer  monitors  the  flow  rate  for 
Him  cooling. 

A  photograph  of  the  chamber,  the  row  of  film  cooling  holes  and 
vortex  generator  is  shown  in  Figure  3.b.  The  chamber  is  constructed  of  1.27 
cm  plexiglass,  with  outside  dimensions  of  0.305  x  0.508  x  0  457  m  The 
internal  structure  consists  of  three  thin  metal  plates  0.381  x  0.508  m. 
starting  5  08  cm  from  the  bottom  and  procee  ling  up  at  5.08  cm  intervals 
Two  silicon  rubber  heaters,  0.381  x  0.483  m,  120  volts,  are  separately  placed 
over  the  bottom  two  plates.  The  heaters  are  controlled  through  a  Powerstat 
variable  autotransformer,  type  136.  The  top  surface  contains  13  plexiglass 
injection  tubes  8  cm  Ions  each,  with  an  inside  diameter  of  0.952  cm  ( 3/8  in  i 
with  a  1/d  ratio  of  8.42.  The  13  injection  nozzles  are  incline  at  angle  of  30°  , 
with  a  three  diameter  spanwise  spacing  between  center  lines.  The  middle 
tube  is  located  on  the  center  line  of  the  test  surface. 

Three  pressure  taps,  positioned  at  the  center  of  the  front  and  two 
side  faces  of  the  injection  plenum  chamber,  are  used  to  measure  Poc  -  P^ 
Three  0  254  mm  diameter  copper-constantan  wire  thermocouples  with 
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welded  joints  are  placed  at  different  locations  inside  the  plenum  chamber 
measure  Top 

2.  Qualification  and  Performance. 

Extensive  qualification  of  the  inaction  system  was  conducted  bv 
|oseph,  [Ref.  3:  pp.  23-26],  Results  show  that  the  uniformity  of  the  olenum 
chamber  pressure,  Poc,  was  satisfactory  over  the  range  of  the  injection  con¬ 
ditions  with  differences  of  about  1%  in  the  spanwise  direction  and  a  ma  o- 
mum  of  4%  occurring  for  only  one  case  at  low  flow  rate  of  0.327E-4  mVs 
This  is  equivalent  of  a  blowing  ratio  M  of  0  004  at  lioo  =  10  rn/s 

The  plenum  produces  a  reservoir  of  air  at  an  elevated  temperature 
and  pressure,  which  is  near  stagnation  conditions.  The  temperature  at  the 
exits  of  the  injection  tubes,  Trc.  is  different  from  the  temperature  of  the 
plenum  of  the  injection  chamber,  Top,  due  to  conduction  through  the  tubes 
surlace  to  the  surrounding  air.  It  is  more  convenient  to  measure  fop, 
whereas  Trc  is  needed  to  calculate  the  injection  parameters  Thus  a 
relation  between  Trc  and  Top  was  needed  and  found  from  experiment  to  be 
T0c  =  l.'i55  ToPIJS6S  Because  of  the  low  velocities  employed,  land  negligible 
viscous  dissipation).  Toc  -  Trc  within  a  fraction  of  a  degree. 

In  order  to  determine  injection  parameters,  the  following  quantities 
must  be  measured,  Poo,  Too  Vc,  Poc,  T0c.  and  A,  the  area  normal  to  the  flowr 
of  the  injection  holes,  designed  to  be  9.2633E-4  m2  The  coolant  velocity  is 
then  given  by: 
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The  static  density  is  estimated  using  : 


Pc  = 


RTC 


>eqn.  l.l> 


The  mass  flux,  mc,  is  the  product  of  Uc  and  pc  .  To  calculate  the  isentropic 


nass  flow.  Pei  and  Uc;  are  found  using 


Pci  = 


P,o 


RTor 


' eqn.  2. 


and 


"v 


2(Poc  "  Poo) 


Pci 


(eqn.  2.4 ) 


The  product  of  pCi  and  Uci  for  compressible  flows  may  be  given  as  : 


0  5 


(rcuc)rpJ  if 


RTocll- 


(eqn.  2.5) 


Discharge  coefficients,  C<|,  are  then  estimated  using 


Cd  = 


PcUc 

(pcUc)i 


i.eqn.  2,6' 


In  order  to  check  the  injection  system,  discharge  coefficients  were 
measured  as  volumetric  flow  rale  was  changed.  Results  are  given  in  Figure 
4  as  a  function  of  Reynolds  number  for  13  injection  hole  locations  and  7 
injection  hole  locations.  Discharge  coefficients  range  between  0.5  and  0.730. 


and  increase  with  Revnolds  number.  For  Re  3  x  UP,  discharge  coefficients 
collapse  closelv  around  a  value  of  approximately  0  72  Because  these  re¬ 
sults  are  consistent  with  those  of  [Ref.  8],  satisfactory  injection  system 
performance  is  indicated. 

C.  HKAT  TRANSFHK  SURFACE. 

The  heat  transfer  surface  was  designed  and  developed  to  provide  a  con¬ 
stant  heat  flux  over  its  area.  The  average  surface  temperature  may  be 
adjusted  and  maintained  from  ambient  up  to  about  60  2C  The  plate  was 
constructed  so  that  its  upward  facing  part  is  adjacent  to  the  wind  tunnel  air 
stream,  with  minimal  heat  loss  by  conduction  from  the  sides  and  beneath  the 
test  surface.  The  plate  itself  has  been  instrumented  to  measure  tempera¬ 
tures  with  thermocouples  placed  just  beneath  the  foil  surface.  A  film  of 
liquid  crystals  is  sprayed  over  a  layer  of  black  paint  painted  on  the  foil. 

1.  Description, 

A  photograph  of  the  heat  transfer  surface  is  shown  in  Figure  5  a 
The  design  is  based  on  ones  used  at  the  University  of  Minnesota  (Ref  13! 
and  [Ref.  14],  and  provides  more  uniform  heat  flux  and  better  spatial  resolu¬ 
tion  of  temperature  than  the  surface  used  by  Joseph,  [Ref.  3:  pp.  27-36],  It 
consists  of  a  thin  stainless  steel  foil,  AISI  302  full  hard,  0.2032  mm  x  1.3  m  x 
0  467  m,  painted  flat  black  with  7  layers  of  liquid  crystals.  Attached  to  the 
under  side  of  the  surface  are  126  copper-constantan  thermocouples.  120  of 
these  have  flattened  tips  and  were  manufactured  by  Marchi  associates  ftyne 
MA  3%T-30-96-FEP-FJ),  the  remaining  six  have  round  junctions  and  were 
manufactured  by  Omega  Engineering,  Inc.  Thermocouple  lead  wires  are 


located  in  grooves  cut  into  a  triple  sheet  ol'  0.254  mm  ( 10  mil)  thick  double 
sided  tape,  manufactured  by  the  3M  Company),  as  shown  in  Figure  5  b.  The 
grooves  are  then  filled  with  RTV.  A  thin  foil  constant  heat  flux  heater  1 .0 
mm  x  1.143  m  x  0.457  m,  120  V/ 1 500  W,  manufactured  by  Marchi 
Associates,  is  attached  to  the  tape  with  Electrobond  epoxy  Beneath  the 
heater  is  a  12.7  mm  !  1/2  ini  thick  lexan  sheet,  followed  by  a  25  4  mm  oi 
foam  insulation,  82.55  mm  thick  styrofoam,  three  sheets  of  254  mm  eacn 
lexan  and  one  sheet  of  9.53  mm  thick  balsa  wood,  as  shown  in  Figure  7. 

Around  the  edges  of  the  foil,  grease  was  inserted  between  it  and  me 
plexiglass  frame  to  fill  any  small  gaps  resulting  from  thermal  expansion.  In 
addition,  both  the  upstream  edge  and  the  trailing  edge  were  taped  to  the 
bottom  wall  of  the  test  section.  Since  additional  vertical  movement  of  the 
foil  above  the  bottom  wall  of  the  test  section  occurs  due  to  thermal  expan¬ 
sion  during  heating,  the  level  of  the  surface  is  adjustable  and  maintained 
level  with  the  test  surface  by  adjusting  screws  in  the  plexiglass  frame  sup¬ 
porting  the  heat  transfer  surface  from  below.  During  heat  transfer  tests,  the 
top  surface  of  the  foil  remained  remarkably  flat  and  smooth  with  minimal 
surface  irregularities.  The  surface  temperature  is  controlled  by  adjusting 
input  voltage  to  the  heater  using  a  Standard  Electrical  Product  Co.  variac, 
type  3000B.  With  this  type  of  heat  transfer  surface,  good  resolution  of  tem¬ 
perature  was  achieved  without  hot  spots. 

Thermocouples  are  placed  on  the  surface  as  shown  in  Figure  6.  In 
each  of  the  six  rows,  21  thermocouples  are  located  1.27  cm  apart  from  each 
other  The  first  row  of  thermocouples  is  located  at  5  cm  from  the  leading 
edge  of  the  test  surface,  the  second  is  10  cm  further  downstream  with 
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respect  to  the  first,  the  third  is  15  cm  apart  downstream  from  the  second 
row,  and  the  rest  of  the  thermocouple  rows  are  20  cm  apart  from  each  other 


The  present  test  surface  was  the  second  one  constructed  lor  this 
experiment.  The  first  one  could  not  be  used  since  it  had  hot  spots  located 
under  the  thermocouple  lead  wires  where  the  surrounding  liner  was  not 
property  attach  to  the  heater.  In  the  second  surface,  this  problem  was 
remedied  by  increasing  the  thickness  of  the  liner 

To  test  the  second  heat  transfer  surface,  a  Huges  Probeye  Thermal 
Series  4000  Video  System  ,  consisting  of  a  infrared  and  video  camera  with 
display  screen,  was  used.  Surface  temperatures  were  measured  as  the  test 
surface  was  heated  outside  the  wind  tunnel,  since  the  system  could  not  see 
through  wind  tunnel  walls.  With  this  system  temperature  variations  as 
small  as  1  QC,  can  be  measured.  The  surface  was  observed  under  three  ope¬ 
rating  conditions  :  natural  convection  with  a  surface  temperature  of  ap¬ 
proximately  35  QC,  forced  convection  from  the  leading  edge  with  a  surface 
temperature  of  about  40  QC,  and  force  convection  from  the  trailing  edge, 
with  a  surface  temperature  of  about  40  2C.  Results  showed  that  most  of  the 
heat  transfer  surface  temperatures  were  spanwise  uniform  within  a  fraction 
of  a  degree  for  all  three  tests.  The  only  exceptions  were  several  small  cool 
spots  with  temperatures  about  l  QC  lower  than  the  rest  of  the  plate,  located 
between  rows  of  thermocouple  lead  wires  near  the  edges  of  the  plate. 

A  second  test  using  the  liquid  crystals  was  performed  to  further 
qualify  the  test  surface.  The  liquid  crystals,  manufactured  bv  Davis  Liquid 
Crystals,  Inc.,  are  rated  from  30  to  35  2C.  Temperature  variations  as  small  as 
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1.2  aC  can  be  measured  with  an  uncertainty  of  t  2  yC  The  test  was  con¬ 
ducted  with  natural  convection  with  a  surface  temperature  of  approximated 
32  aC.  Results  were  consistent  with  the  infrared  test  with  temperatures 
variations  not  greater  than  1.2  across  most  of  the  surface.  As  for  the 
early  test,  small  hot  spots  were  evident  near  some  thermocouple  lead  wire 
locations. 

3 .  Energy  Balance. 

An  energy  balance  was  performed  to  determine  the  heat  loss  by 
conduction  from  the  heat  transfer  test  surface  used  to  obtain  final  results 
During  the  energy  balance,  heat  loss  by  radiation  and  convection  were  pre¬ 
vented  since  the  metal  foil  surface,  ordinarily  exposed  to  convection  in  the 
wind  tunnel,  was  covered  with  three  layers  of  25.4  mm  thick  foam  insula¬ 
tion.  For  the  energy  balance,  and  for  all  wind  tunnel  testing,  foam  insulation 
was  also  placed  around  the  sides  of  the  test  surface  located  below  the  wind 
tunnel  convection  surface.  To  estimate  heat  loss  through  the  insulation 
placed  on  top  of  the  foil  surface,  program  ENERGB  was  employed  with  an 
algorithm  involving  the  one  dimensional,  linear  form  of  Fourier  s  conduction 
equation: 

AX 

qw  =  1CA—  (eqn.  2.7) 

For  the  insulation,  1C  is  0.4  W/m  aC,  A  is  0.4897  m2,  AX  is  0.0254  ra,  and  AT 
is  the  temperature  drop  in  the  X  direction  in  aC.  Heat  conduction  through  the 
bottom  and  sides  of  the  heat  transfer  device  is  then  given  by: 
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qc  =  VI  -  qv 


(eqn.  2.8) 
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Here.  VI  is  the  power  mio  the  test  piate,  and  qw  is  the  conduction  loss 

through  the  top  insulation.  Tests  were  made  at  live  power  levels  1 1.02.  ; 

13.0,  15.25,  16.52  and  18.2  Watts,  chosen  to  give  conduction  losses  at  the  \ 

f 

same  levels  as  experienced  under  normal  operating  conditions.  Figure  8  • 

snows  qc  vs.  Iv  -  T*mb  results,  where  is  the  average  piate  surface  < 

temperature.  A  second  order  polynomial  was  fitted  to  this  data  in  order  to  ; 

predict  conduction  losses  during  heat  transfer  measurements: 

( 

t 

qc  =  0.683  +  0.954Tdiff  -  0.01 6Tdiff2  ieqn.  2.9? 

where  Tdiff  =  Tw  -  Tamb-  This  equation  is  valid  over  a  range  of  Tw  -  Tamb 
from  10  to  30  QC.  When  exposed  to  convection  in  the  wind  tunnel,  conduc¬ 
tion  losses  are  only  1.5  to  2.5%  of  the  total  power,  and  thus,  a  25%  error  in 
the  estimate  of  conduction  tosses  will  cause  less  than  a  0.5%  error  in  the 
estimate  of  the  heat  transfer  by  convection. 

Radiation  losses  were  estimated  using  two  different  approaches. 

For  the  first: 
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(eqn.  2.10) 
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drad  =  Xdij 


(eqn.  2.11) 
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[Ret.  1 6 1.  The  view  factors,  Fij,  for  the  top  and  each  of  the  side  wails  were 
estimated  to  be  0.54  and  0.23.  respectively  [Ref.  16],  where 
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■  eqn.  2.13 


From  these  two  methods,  the  radiation  losses  for  an  average  plate  tempera¬ 
ture  of  40  QC  were  estimated  to  be  55  Watts,  approximately  8. 5*  of  the  total 
power  into  the  heat  transfer  surface. 

4.  Contact  Resistance  Temperature  Pro 


Local  surface  temperatures  were  measured  using  thermocouples 
placed  in  contact  with  the  underside  of  the  metal  foil  surface.  The  junctions 
of  these  thermocouples  were  held  in  place  next  to  the  foil  with  three  lavers 
of  0.254  mm  thick  double-sided  lining  tape  and  RTV  silicon  rubbei  epoxy. 
During  the  heat  transfer  tests,  temperatures  measured  by  the  thermocouples 
were  greater  than  those  of  the  surface  of  the  plate,  due  to  thermal  contact 
resistance  between  the  thermocouples  and  the  foil  and  conduction  through 
the  foil.  The  resulting  temperature  difference  AT  is  given  by: 


AT  -  qj 


!.  hc  A  LA 


i  eqn.  2.11 


where  AX/XA  accounts  for  conduction  across  the  foil  thickness  and  l/hcA 
accounts  for  the  thermal  contact  resistance,  q  is  the  heat  flux  through  the 
foil.  The  contact  resistance  is  highly  dependent  on  the  contact  pressure  as 
well  as  the  area  of  the  surfaces  in  contact. 

The  value  of  contact  resistance  used  in  the  present  studv  was  the  same 
as  used  Oy  Joseph  [Ref.  51.  in  h is  study  liquid  crystals  were  applied  to  the 
surface  to  measure  its  temperature  distribution  during  convection  tests,  so 
that  AT  in  eqn.  2.13  could  be  determined.  Stanton  numbers  from  baseline 
tests  matched  expected  correlations  after  accounting  for  contact  resistance. 
The  empirical  relationship  for  turbulent  boundary  layers  at  constant  free 
stream  velocity,  along  a  flat  plate  with  constant  wail  heat  flux  and  unheated 
starting  length  of  £  =  1.10  m.  (Ref.  Ml  is  given  by 
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Data  in  Figures  10  and  13  show  that  this  equation  matches  the  data  well  for 
5  x  105-  Re*  2.0  x  106  From  these  measurements,  the  two  terms  in 
brackets  in  equation  2.  13  were  estimated  to  be  0.014  2K/Watt.  The  same 
value  was  used  for  all  thermocouples.  However,  on  the  present  test  plate 
contact  pressure  and  the  area  of  the  thermocouples  in  contact  with  the  sur¬ 
face  varied  slightly  from  one  thermocouple  to  another.  Because  contact 
resistance  for  each  individual  thermocouple  may  vary  from  this  value,  small 
deviations  in  the  spanwise  heat  transfer  coefficients  and  Stanton  number? 
result,  which  were  independent  of  flow  conditions.  In  the  present  study,  the 
effects  of  these  small  variations  were  minimized  by  presenting  results  for 
local  conditions  in  terms  of  Stanton  number  ratios . 


D.  TEMPERATURE  MEASUREMENT. 


All  thermocouples  employed  in  the  present  study  are  type-T,  copper- 
constantan  thermocouples.  120  of  these  were  manufactured  by  Marchi 
Associates.  Inc  and  attached  beneath  of  the  heat  transfer  surface.  These 
thermocouples,  have  a  flattened  junctions,  so  that  they  may  be  placed  in 
good  contact  with  the  surface.  One  of  these  thermocouples  was  calibrated 
using  a  temperature  regulated  bath  consisting  of  liquid  nitrogen  and  electric 
heaters  and  a  platinum  resistance  temperature  reference  (±  0.01  &C).  The 
calibration  was  performed  over  the  temperature  r.nge  of  0  -  45  32  °C.  A 
second  order  was  used  to  convert  microvolts  to  temperature  : 

T  -  0.018205  +  0.025846E  -  0.00000058 IE2  (eqn.  2.16) 

where: 

E-  microvolts  x  106  Volts  (eqn.  2.17) 

This  calibration  was  used  for  all  Marchi  Associates  thermocouples  attached 
to  the  plate,  since  all  thermocouples  indicated  very  similar  behavior  as 
temperatures  were  changed. 

Thermocouples  manufactured  by  Omega  Engineering  Co.  were  used  to 
measure  six  plate  temperatures,  as  well  as  temperatures  of  the  plenum  air, 
freestream  air,  and  the  boundary  layer  as  traverses  were  made.  The  calibra¬ 
tion  of  Joseph,  (Ref.  3.:  pp.  37],  was  first  verified,  and  then  used  for  all  of 
these  measurements  except  the  traverses: 


T  =  26.573E  -  1.937E2  +  0.998E?  -  0  26 1 


'eqn  2.18) 


Here, 

E  =  millivolts  x  1000  (eqn.  2.19) 

Temperatures  surveys  of  T  -  Tco,  were  performed  using  Hie  automated 
traversing  device.  For  this,  two  thermocouples  manufactured  by  Omega 
Engineering  were  used,  after  calibration  using  the  technique  mentioned 
above.  Simultaneous  measurement  of  free  stream  temperature  was  made 
for  every  boundary  layer  location  to  minimize  scatter  due  to  wind  tunnel 
temperature  drift.  For  the  thermocouple  mounted  in  the  probe  that  travels 
the  flow  field,  the  following  equation  was  used  to  convert  voltage  to 
temperature: 

T  -  0.1836  +  0.025667E  -  0.0000004882E2  (eqn.  2.20) 

For  the  thermocouple  used  to  measure  the  free  stream,  the  equation 
employed  was  : 

Too  =  0.07832  +0.023 1054E  -  0.0000006786E2  (eqn.  2.21 ) 

As  before,  E  -  microvolts  x  106.  The  automated  traversing  mechanism  has 
two  degrees  of  freedom  which  allowed  measurement  of  a  plane  the  flow 
field.  Each  survey  consisted  of  800  probe  locations,  covering  an  area  of  12 
cm  x  22  cm,  Both,  spanwise  and  vertical  traversing  blocks  are  mounted  on  a 
20-thread  per  inch  drive  screw  and  two  ground  steel  case  hardened  steel 
guide/support  shafts.  Each  drive  shaft  is  directly  coupled  to  a  SLO-SYN  type 


M092-FD310  steeping  motor.  Motors  are  controlled  by  a  MIT  AS  two-axis 
Motion  Controller.  Both  the  motors  and  the  controller  are  manufactured  by 
The  Superior  Electric  Co.  The  MIT  AS  controller  comes  equipped  with  2k 
bytes  of  memory  and  a  MC68000,  1 6 -bit  microprocessor  which  allows  the 
user  to  program  the  start,  stop,  duration,  speed. acceleration  and  deceleration 
of  the  steeping  motors. 

E.  DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system  was  designed  to  rapidly  measure  thermo¬ 
couple  voltages  and  convert  them  to  temperature  in  degrees  C.  Using  these 
temperatures  along  with  user  supplied  information  on  ambient  conditions, 
freestream  conditions,  power  input,  and  flow  rates  into  the  injection 
chamber,  the  system  calculates  free  stream  velocity,  densitv,  local  heat 
transfer  coefficients,  Stanton  numbers,  spanwise  averaged  Stanton  numbers 
and  injection  parameters  such  as  blowing  ratio  and  discharge  coefficient. 

1.  Hardware 

A  Hewlett-Packard  Series  300,  Model  9836S  computer,  equipped 
with  a  MC68000,  8  MHz  16/32-bit  processor,  dual  5-1/4  inch  floppy  disk¬ 
drives,  and  1M  bytes  of  memory  RAM,  was  used  in  the  study.  The  computer 
was  used  to  process  signals  from  the  data  acquisition  system  as  well  as  the 
information  supplied  by  the  user  in  order  to  process,  store,  display  and  print 
results.  An  HP  Think  Jet  printer  was  used  to  print  data,  and  an  HP  7470  two 
pen  plotter  was  employed  for  graphics  Voltages  from  the  thermocouples 
are  read  by  an  HP-3497A  Data  Acquisition/Control  Unit  with  an  HP-34^8 A 
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Extender.  The  unit  communicates  with  the  computer  through  a  HP-X2(i7';~  \ 

!  nten  ace 

2.  Software. 

Programs  STANFC1,  STANFC2  and  ACQTPKO  were  developed  to 
process  temperature  and  Stanton  number  data,  program  ENERGB  was  deve- 
ioped  to  estimate  the  conduction  losses,  and  programs  PTSLC  PiSi'A  , 
PLSTRVOR.  PLSTRTIO.  PLSTRFC.  PLSTVV,  SL'RFCONT  and  PLOTRl'N  were 
developed  tor  plotting  data.  All  these  programs  are  listed  in  Appendix  C. 
The  data  files  created  by  the  heat  transfer  programs  are  listed  in  Appendix 
D. 

Programs  STANFC1  and  STANFC2  are  modified  versions  of  STDAT1. 
which  was  developed  by  Ligrani,  Ortiz  and  Joseph.  [Ref  3.:  PP  41]  Program 
ST  A  NFC  1  prompts  the  user  if  film  cooling  is  being  used.  If  the  answer  is 
affirmative,  the  program  prompts  the  user  to  enter  the  percentage  of  flow  to 
the  injection  chamber  from  the  rotometer  and  the  pressure  difference 
between  the  injection  chamber  and  the  static  pressure  in  the  wind  tunnel 
These  inputs  are  transformed  to  SI  units  and  stored  in  a  data  file  named 
F I  LOT".  The  program  continues  by  prompting  the  user  for  the  stagnation 
pressure  of  the  free  stream  (in  inches  of  H2O),  the  ambient  pressure  <  in  Hg 
the  current  (Amps)  and  voltage  supplied  to  the  heater.  The  program  then 
reads  the  thermocouple  voltages  and  converts  them  into  temperatures 
These  are  then  stored  in  a  data  file  named  "TDAFC".  After  all  temperatures 
have  been  calculated,  the  free  stream  density  and  velocity  are  calculated 
Parameters  previously  calculated  such  as  the  free  stream  density,  tree 
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stream  velocity  ambient  pressure,  and  the  power  in  are  stored  in  a  data  tile 
named  ID  AFC 

Program  STANFC2,  processes  the  data  tiles  TDArC,  1DAFC  and  FILDT 
The  program  begins  by  prompting  the  user  whether  a  vortex  generator  is 
being  used  or  not.  then  it  accounts  for  conduction,  radiation  losses  and 
contact  resistance  It  continues  bv  calculating  local  heat  transfer  coet haems 
and  Stanton  numbers,  average  spanwise  Stanton  numbers  and  Reynolds 
number  based  on  downstream  distance  These  data  are  printed  out  and  two 


data  files,  STRFCY  and  STAY'.  are  created.  In  STRFC"  the  local  heat 
transfer  coefficients  and  Stanton  numbers  are  stored,  along  with  their 
spanwise  position  and  downstream  positions  along  the  test  plate.  Spanwise 
averaged  Stanton  numbers  with  their  corresponding  down  stream  Reynolds 
numbers  are  located  in  STAY  1'he  last  section  of  the  program  includes  a 
subroutine  to  calculate  film  cooling  parameters,  such  as  discharge  coeffi¬ 
cients,  density  ratios,  mass  flux  ratios,  momentum  flux  ratios  and  blowing 
rate 

Program  ACQTPRO.  developed  by  Ortiz,  acquires  temperatures  from  a 
thermocouple  probe  mounted  on  the  automatic  traversing  device  and  from 
another  thermocouple  that  senses  the  freestream  temperature.  This  pro¬ 
gram  was  used  to  create  the  data  to  plot  the  temperature  surveys  of  T  -  Too. 
The  program  begins  by  prompting  the  user  for  the  downstream  distance 
from  the  boundary  layer  trip,  the  number  of  points  in  the  spanwise  direction 
where  temperatures  are  going  to  be  measured,  the  number  of  points  tn  the 
'vertical  direction  the  spanwise  resolution  in  mco  is,  the  vertical  resolution  in 
inches  and  the  initial  position  of  the  probe  with  respect  to  the  plate  in  both  /. 
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and  Y  coordinates,  in  inches.  A  matrix  of  data  points  is  then  computed 
Preesiream  temperature  is  calculated  and  ambient  conditions  are  input  Tv 
program  then  enters  a  loop  which  samples  each  temperature  25  times  per 
probe  position  and  5  times  the  free  stream  temperature.  Upon  acquiring 
these  temperatures  and  obtaining  their  respective  average,  T  -  To  is  calcu¬ 
lated  for  each  probe  location.  The  temperature  difference  along  with  their 
respective  coordinates  are  stored  in  a  matrix.  At  the  end  of  the  data 
collection  run.  these  values  are  read  into  a  data  lile  named  TPKO  on  a 
floppy  disk  These  data  are  plotted  using  program  TPKOPUN 


Ill  EXPERIMENTAL  RESULTS 


A.  BASE  LINE  MEASUREMENTS 

Heat  transfer  measurements  were  made  at  free  stream  velocities  of  10 
m/s  and  15  m/s,  without  vortex  and  without  film  cooling.  These  were  done 
in  order  to  determine  how'  heat  transfer  from  the  plate  compares  with 
existing  correlations.  These  measurements  validated  and  qualified  the  per¬ 
formance  of  the  heat  transfer  plate  and  measurements  procedures  used 
Results  are  given  in  Figures  9-11  for  10  m/s,  and  in  Figures  12-13  for  15 
m/s. 

Figure  9  shows  the  spanwise  uniformity  of  the  local  heat  transfer  coeffi¬ 
cient  of  the  test  surface  at  10  m/s.  Except  for  row  1,  the  spanwise  unifor¬ 
mity  is  very  good  with  maximum  variation  of  10%  (based  on  the  average  for 
a  given  row).  These  small  variations  are  probably  due  to  slight  differences 
ir.  the  spatial  uniformity  of  the  heat  transfer  test  surface,  especially  diffe¬ 
rences  in  the  conduction  contact  resistance  between  different  thermocouples 
and  the  metal  foil,  which  comprises  the  top  of  the  heat  transfer  plate  Larger 
spanwise  variations  for  row  1  are  believed  to  be  due  to  multi-dimensional 
heat  transfer  by  conduction  through  the  leading  edge  and  through  the  front 
corners  of  the  test  plate  in  addition  to  the  contact  resistance.  These  varia¬ 
tions  are  not  as  large  as  those  observed  by  Joseph  [Ref.  3:  pp  45,  77-78) 

The  spanwise-averaged  Stanton  Numbers  for  10  m/s  are  shown  in  Fi¬ 
gure  10.  Six  data  points  are  shown,  where  each  corresponds  to  a  different 
thermocouple  row  The  data  shows  excellent  agreement  with  the  empirical 


equation  for  turbulent  boundary  layers  on  a  flat  plate  at  constant  tree 
stream  velocity,  with  constant  heat  flux  and  unheated  starting  length  of  ■  - 
1.10  m  [Ref.  151. 


A  survey  of  T  -  TA=  in  degrees  Centigrade  was  obtained  at  X  =  1.  m. 
using  a  calibrated  copper-constantan  thermocouple,  positioned  using  the 
automated  traversing  device.  Results  are  shown  in  Figure  11.  for  LAo-10 
m/s  with  no  vortex  and  no  film  cooling.  In  this  figure,  results  for  an  area 
normal  to  the  flow  of  1 1  cm  x  20  cm.  are  given,  as  measured  using  800 
probe  positions.  The  survey  evidences  good  spanwise  uniformity  of  the 
temperature  boundary  layer. 

Figure  12  and  13  show  the  spanwise  variations  of  the  heat  transfer 
coefficient  and  averaged  Stanton  numbers  as  function  of  Reynolds  number 
for  15  m/s  free  stream  velocity.  Results  are  consistent  with  those  obtained 
for  10  m/s.  In  Figure  12,  the  data  show  good  spanwise  uniformity,  except 
for  the  first  row,  where  variations  are  slightly  larger  than  for  the  10  m/s 
case.  In  Figure  13,  the  data  again  show  excellent  agreement  with  equation 
(3.1). 

B.  BOUNDARY  LAYER  WITH  SINGLE  VORTEX 

Heat  transfer  measurements  were  made  at  free  stream  velocities  of  10 
m/s  and  15  m/s  with  an  embedded  longitudinal  vortex  and  no  film  cooling. 
These  were  done  to  further  qualify  the  heat  transfer  surface  for  spatially  re- 


solved  heat  transfer  measurements.  To  produce  the  vortex,  vortex  generator 
*2,1  see  Joseph  [Ref.  3:  PP  76j),  was  positioned  0.479  m  downstream  from  me 
boundary  layer  trip,  and  4.79  cm  in  the  positive  Z-direction,  as  shown  in 
Figure  2. 

The  spanwise  variation  of  local  St/Sto  for  10  m/s  and  15  m/s  are  pre¬ 
sented  in  Figures  14  and  15,  respectively.  Here,  Sto  is  the  local  Stanton  num¬ 
ber  without  an  embedded  vortex.  In  both  of  these,  the  effect  of  the  vortex 
on  wall  heat  transfer  is  evident.  The  Stanton  number  ratios  approach  l 
away  from  the  vortex.  As  Z  decreases  from  W5cm,  the  ratios  increase  until 
maximum  values  of  1.2  and  1.25  are  reached  for  10  and  15  m/s,  respec¬ 
tively.  These  maxima  correspond  to  locations  on  downwash  sides  of  vor¬ 
tices,  where  secondary  flows  result  in  boundary  layers  which  are  locally 
thinner  than  at  other  locations  For  -5.0<Z<.0.0  cm  a  large  gradient  ot  heat 
transfer  exists  at  each  downstream  location.  The  location  of  this  gradient 
moves  to  smaller  Z  as  the  vortex  develops  downstream.  For  smaller  Z, 
Stanton  number  ratios  for  10  and  15  m/s,  reach  minima  of  0.92  and  0.90  at 
locations  which  correspond  to  the  upwash  sides  of  vortices.  Here,  the 
boundary  layer  is  locally  thicker  than  at  other  locations. 

The  results  in  Figures  14  and  15  are  in  excellent  quantitative  agree¬ 
ment  with  those  of  Joseph  [Ref.  3:  pp.  87-88],  for  the  same  experimental 
conditions,  (The  sign  of  the  Z  coordinate  in  this  experiment  is  opposite  to 
the  one  employed  by  Joseph).  The  results  in  14  and  15  also  show  qualita¬ 
tive  agreement  with  the  measurements  of  Eibeck  and  Eaton  [Ref  6],  where 
small  quantitative  differences  are  a  result  of  different  vortex  generator 
geometries. 
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C.  BOUNDARY  LAYER  WITH  FILM  COOLING 


In  order  to  further  qualify  the  heat  transfer  surface  and  the  film  cooling 
injection  system,  measurements  were  made  with  film  cooling  without  em¬ 
bedded  vortex.  Results  were  obtained  at  10  m/s  ,  at  9  =  1.614.  1.609  and 
;  y.)  with  blowing  ratios  of  0  68.  0  08  and  1  26  respectively  and  at  !  8  m/- 
at  9  =  1.393  and  1.886  with  blowing  ratios  of  0.47  and  0.86  respectively 
Results  are  presented  in  Figures  16-30. 

Figures  16-18  show  the  spanwise  variations  of  St/Sto  at  10  m/s,  for 
blowing  ratios  of  0.68,  0.98  and  1.26,  respectively.  The  first  data  set  was 
obtained  using  13  injection  holes,  which  produce  a  film  wide  enough  (23.4 
cm)  to  cover  the  entire  span  of  the  measuring  heat  transfer  surface  Because 
of  the  limitations  on  the  secondary  air  supply,  9  and  7  injection  holes  were 
required  to  produce  blowing  ratios  of  0.98  and  1.26,  respectively  In  these 
cases  the  coolant  was  injected  such  that  it  covered  a  portion  of  the  center 
span  of  the  heat  transfer  test  surface.  Consequently,  Stanton  number  ratios 
are  higher  near  the  edges  of  the  plate,  as  shown  in  Figures  17  and  18  The 
latter  Figure  for  only  7  injection  holes,  shows  this  effect  to  be  particularly 
significant,  however,  in  spite  of  this  a  large  spanwise  uniform  area  exists  in 
the  center  portion  of  the  plate  where  adequate  spanwise  averages  may  be 
obtained.  Apart  from  this,  the  spanwise  uniformity  of  Stanton  number  ratios 
in  Figures  16-18  is  excellent  with  a  maximum  deviation  of  about  10  percent. 

Figures  19  and  20  show  the  spanwise  variations  of  the  Stanton  number 
ratios  at  13  m/s.  Results  for  blowing  ratios  of  0.47  and  0  86  were  obtained 
using  1  j  and  7  injection  holes  respectively.  Spanwise  variations  in  Figures 


19  and  20  are  larger  than  in  figures  lb- 18.  though  qualitative  trends  are 
verv  similar. 


A  summary  ol'  spanwise  averaged  St/Sto  as  dependent  upon  Reynolds 
number  is  given  in  Figure  21.  Graphs  showing  individual  curves  are  pre¬ 
sented  in  Figures  22-  26,  where  averaged  St/Sto  are  given  for  blowing 
ratios  of  0.b8.  0.98  and  1.26  at  10  m/s.  and  0.47  and  0.86  at  IS  m/s  respec- 
tivelv.  Referring  to  Figure  21,  the  lowest  St/Sto  are  observed  for  \l  -  o  ai 
10  m/s  and  for  M  =  0.47  at  15  m/s.  For  each  individual  data  set  the  lowest 
ratio  is  present  for  the  two  rows  of  thermocouples  nearest  the  film  cooling 
holes,  where  X/d  is  7.35  and  17.48.  Here,  X  is  distance  from  the  down 
stream  edge  of  the  injection  holes,  and  d  the  injection  hole  diameter  The 
data  in  figures  21-26  are  plotted  as  function  of  the  blowing  ratio  in  figures 
27  and  2S.  Results  in  21-28  are  qualitatively  consistent  with  the  fully  tur¬ 
bulent  measurements  of  Goldstein  and  Yoshida,  iRef.  151.  also  for  a  single 
row  of  injection  holes. 

Two  surveys  of  T  -  Too  in  degrees  C  were  obtained  in  the  film  cooled 
turbulent  boundary  layer  at  X  =  1.48  m,  X/d  =  41.94.  Both  survevs  were 
made  at  a  free  oircam  velocity  of  10  m/s,  with  M  =  0.98  and  a  coolant  infec¬ 
tion  temperature  of  5 1  QC.  In  one  case,  the  plate  was  heated  to  40  OC  (Figure 
30  ),  and  in  the  other,  the  plate  was  maintained  at  the  free  stream  tempera¬ 
ture  (Figure  29  ).  The  latter  experimental  arrangement  with  unheated  plate 
allowed  the  determination  of  the  presence  and  distribution  of  heated  film 
coolant  as  indicated  by  higher  temperatures.  In  Figure  30  with  a  heated 
plate,  the  locations  of  coolant  jets  are  evident  as  local  hot  spots  near  the  wall 
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Comparing  these  results  to  ones  in  Figure  1 1.  without  film  cooling  shows  that 
the  thermal  boundary  layer  is  about  twice  as  thick. 


D.  BOUNDARY  LAYER  WITH  SINGLE  VORTEX  AND  FILM  COOLING 

In  this  part  of  the  study  the  effect  of  an  embedded  vortex  on  a  film 
cooled  turbulent  boundary  layer  was  studied.  This  section  of  the  experiment 
was  conducted  in  four  different  parts.  For  all  the  parts  of  the  experiment 
the  same  vortex  generator  was  employed.  The  first  consisted  in  repeating 
Joseph  data.  (Ref.  3.  :  pp.  97-104):  Uoo-  10  m/s,  vortex  generator  *2  at  Z  - 
4.79  cm  and  blowing  ratio,  M  *  0.68,  in  order  to  provide  additional  check  on 
measurements,  equipment  and  procedures.  In  the  second  step,  a  blowing 
ratio  of  M  »  0.98  was  employed,  at  U<»  *  10  m/s.  The  position  of  the  vortex 
was  changed  to  three  different  locations  with  respect  to  test  section  center 
line.  Vortex  position  A  corresponds  to  vortex  generator  at  Z  =■  3.52  cm,  vor¬ 
tex  position  B  to  Z  -  4.79  cm,  and  vortex  position  C  to  Z  -  6.06  cm.  In  the 
third  step,  a  freestream  velocity  of  15  m/s,  with  M  -  0.47  were  employed, 
in  order  to  study  the  effect  of  blowing  ratio.  In  the  fourth  part,  a  blowing 
ratio  of  M  -  1.26,  at  Uoo  =  10  m/s  was  employed,  in  order  to  study  the  effect 
of  high  blowing  ratio.  Here,  information  was  also  obtained  on  changes  in 
boundary  layer  structure  which  occurred  for  higher  blowing  ratios. 


Joseph  data. 

These  heat  transfer  measurements  were  obtained  to  compare  with 
those  of  Joseph,  [Ref.  31  He  first  showed  that  embedded  vortices  cause  sig¬ 
nificant  changes  in  heat  transfer  in  film  cooled  turbulent  boundary  layers. 
This  finding  is  also  evident  from  the  present  study.  To  obtain  the  present 


data  set,  vortex  generator  *2  was  placed  at  position  B,  and  all  15  miecuon 
holes  were  used  for  him  cooling.  Results  are  presented  in  Figures  3  la. -51  i 
and  in  Figure  32.  Overall  effects  are  now  discussed.  The  embedded  vortex 
produces  a  thicker  boundary  layer  in  the  upwash  side  thus  augmenting  the 
film  cooling  protection.  Near  the  downwash  side,  the  boundarv  laver  is  io 
cal ly  thinner,  and  the  protection  provided  bv  the  film  cooling  is  minimized 
The  undercooled  region  produced  on  the  upwash  side  vortex  is  verv  ner- 
sistent  not  only  in  the  down  stream  direction,  but  also  in  the  spanwise 
direction.  Here  Stanton  numbers  are  .5  to  5  5  percent  lower  than  those  in 
boundary  layers  with  film  cooling  only 

Figures  3 1  a-f ,  show  the  spanwise  variation  of  the  Stanton  number 
ratios,  St/Sto  with  film  cooling  and  embedded  vortex.  Here,  St  is  the  local 
Stanton  number  and  Sto  is  the  Stanton  number  for  the  same  location  and 
free  stream  velocity  but  without  film  cooling  and  without  embedded  vortex. 
A  St/Sto  =  1  thus  indicates  an  undisturbed  thermal  boundary  layer  As 
shown  in  the  set  of  Figures  31a  -  3  If.  the  lowest  St/Sto  is  0.55  at  X  =  1.15  m. 
X/d  -  7.35,  and  the  maximum  is  1.025  at  X  =  2.00  m.  X/d  =  96.50.  If  these 
results  are  compared  with  those  of  Joseph  [Ref.  3-  :  pp.  97-1031,  a  general 
qualitative  agreement  can  be  observed.  Differences  are  due  to  im¬ 
provements  in  the  present  heat  transfer  test  plate,  especially  better 
spanwise  temperature  resolution  and  more  uniform  heat  flux. 

Figure  32.  shows  surface  contours  of  St/Stf.  Here,  St  is  the  local 
Stanton  number,  and  Stf  is  the  local  Stanton  number  at  the  same  location 
with  film  cooling  but  without  embedded  vortex.  A  very  steep  heat  transfer 
gradient  is  present  which  is  near  the  same  location  as  the  axis  of  the  vortex 
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At  the  upstream  end  of  the  plate  this  gradient  is  located  near  Z  -  -  2  cm  A 
hot  spot  is  present  at  larger  Z.  and  a  region  ot'  high  heat  transfer  is  present  at 
smaller  Z  on  the  upwash  side  of  the  vortex.  These  results  are  in  qualitative 
agreement  with  those  of  Joseph  (Ref.  3.  PP  104],  except  the  most  significant 
hot  spots  are  located  further  downstream  in  the  present  studv 
2  Free  stream  10  m/s.  B low  i  ng_  r  a tio.  M  -  0.98 

In  order  to  achieve  the  blowing  ratio  used  for  this  data.  4  injection 
holes  were  employed  for  film  cooling.  Here,  injection  covers  the  midspan  of 
the  Lest  plate  from  -12  cm  <  Z  s  12  cm.  Vortex  generator  32,  was  placed  hi) 
cm  upstream  of  the  row  of  injection  holes,  at  vortex  position  A,  B  and  C.  to 
investigate  the  effect  of  the  vortex  position  with  respect  to  film  cooling 
injection  holes.  Results  are  presented  in  Figures  33  -  42 
a.  Data  overview. 

Spanwtse  variations  of  St/Sto  for  the  vortex  position  B  are 
shown  in  Figures  33a-33f.  These  data  show  the  same  overall  qualitative 
trends  as  measurements  at  10  m/s  and  M  =  0.68  (Figure  31)  With  embed¬ 
ded  vortex,  the  normalized  Stanton  numbers  in  33a-33f  are  increased  on  the 
downwash  side  and  decreased  on  the  upwash  side  relative  to  the  boundary 
layer  with  film  cooling  only.  The  only  significant  quantitative  difference 
between  these  results  and  those  in  31  are  higher  St/Sto  at  X/d  =  17.48  and 
at  X/d  =  33-59  resulting  from  a  less  effective  film  cooling  at  M  =  0  98. 


b.  Hi'i'eci  01  downstream  development. 


in  Figures  33a.  -  331'.  a  double  peak  oi  St/Sto  is  present  at  X/a 
=  17,48.  and  at  X/d  =  33.59.  As  the  boundary  layer  develops  further  down¬ 
stream,  one  peak  becomes  smaller  (Z  -  -2.0  cm),  and  one  increases  in  mag- 


iarge  peak.  This  St/Sto  peak  increases  in  magnitude  to  become  equal  to 
or  about  34  percent  greater  than  the  St/Sto  with  film  cooling  only.  St/Sii 
data  in  Figure  34.  indicate  that  ,u's  hot  spot  covers  an  area  along  the  center 
and  spanwise  downstream  half  of  the  plate.  Such  behavior  evidences  highly 
three-dimensional  interactions  within  the  boundary  layer. 

Figure  34  shows  the  St/Stf  distribution  for  vortex  position  B 
and  M  -  0.93.  A  steep  heat  transfer  gradient  is  evident  along  the  the  iengtn 
>;f  the  test  surface,  which  corresponds  to  the  path  of  the  vortex  center.  Thu 
path  is  at  an  angle  to  the  streamwise  direction,  moving  to  smaller  Z  as  the 
vortex  convects  downstream.  Compared  to  results  for  M  =  0.68  in  Figure  32. 
the  heat  transfer  gradient  is  steeper,  and  a  larger  region  wThere  St/Stf  is  less 
than  !  is  present  over  the  top  third  of  the  test  plate.  This  follows  since  a 
larger  amount  of  coolant  is  swept  by  the  vortex  at  higher  blowing  ratio. 

For  Z  values  smaller  than  those  along  the  location  of  the  vortex 
center,  Figure  33  shows  values  of  Stanton  ratios  which  are  about  5  per  cent 
lower  than  those  with  film  cooling  only.  The  same  trend  is  evident  in  Figure 
34.  where  almost  1/3  of  the  plate  shows  0.90  -  St/Stf  <  0.98.  These  locations 
are  on  the  upwash  side  of  the  vortex,  where  film  coolant  seems  to  be  pushed 
and  spread  over  a  very  large  portion  of  the  plate.  Similar  phenomena  were 
observed  by  Goldstein  and  Chen.  [Refs.  2  and  71,  in  a  study  of  film  cooling  of 
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a  turbine  blade  with  injection  through  one  and  two  rows  of  holes  in  the 
near -end wall  region 

Figures  35a-  35d  show  the  T  -  Too  temperature  field  as  it  de¬ 
velops  downstream  at  X/d  =  41. 9,  82.9,  109.2  and  147.0.  For  these  tests, 
film  cooling  jets  were  heated  to  51  QC  without  providing  anv  heat  to  the  test 
plate.  Thus  the  temperature  field  shows  how  tluid  from  the  him  cooling 
holes  is  convected  and  distorted  by  the  vortex,  where  higher  temperatures 
indicate  greater  amounts  of  coolant.  The  most  dramatic  effect  evident  from 
Figures  35a-d  occurs  on  the  upwash  side  of  the  vortex,  where  coolant  is 
lifted  away  from  the  thermal  boundary  layer  which  is  ordinarily  about  4  cm 
thick  with  film  cooling.  At  X/d  =  41.9  fluid  affected  by  the  film  cooling  is  6  5 
cm  from  the  surface,  and  at  X/d  =  147.0,  the  temperature  field  is  more  than 
h)  cm  from  the  surface.  Also,  evident  are  cooler  temperatures  with  dow  n¬ 
stream  development  resulting  as  coolant  is  convected  and  diffused.  The 
areas  indicated  by  low  St/Sto  and  St/Stf  values  correspond  to  0.  <  Z  -5  cm, 
where  coolant  seems  lifted  from  the  surface  to  accumulate  in  one  small  area 
Figures  35a- d  clearly  show  how  the  coolant  is  pushed  and  spread  over  an 
area  by  the  vortex,  and  how  this  area  increases  in  size,  as  the  flow  convects 
downstream. 

On  the  downwash  side  of  the  vortex,  secondary  flows  cause 
t'reestream  fluid  to  be  located  very  near  the  wall.  Here,  the  thermal  bound¬ 
ary  layer  is  greatly  thinned  and  very  little  effect  of  the  film  cooling  remains. 

i 

From  Figure  35c  for  X/d  =  82.94,  the  local  hot  spot  is  located  at  to  0  Z  -5 
cm.  which  seems  to  indicate  that  film  coolant  is  moved  to  other  locations  in 
the  boundary  layer  by  the  vortex. 

I 

i 

I 

i 
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In  Figure  36,  the  tempe  -ature  field  is  presented  for  an  expe¬ 
riment  in  which  the  heat  transfer  plate  is  heated  to  40  and  film  coolant  to 
51  °C.  The  temperature  difference  contours  show  the  same  qualitative 
trends  as  in  Figure  35.  In  particular,  the  positions  of  many  temperature 
gradients  and  the  dimensions  of  the  area  affected  by  the  vortex  are  about 
the  same.  Results  in  Figure  36  are  different  from  those  in  35  since  :  (1)  it  is 
not  as  easy  to  discern  which  part  of  the  temperature  field  is  due  to  film 
cooling  only,  (2)  overall  temperature  differences  are  large:,  and  (3)  local  hot 
spots  corresponding  to  film  cooling  locations  are  more  easily  discernabie 
every  3  cm  starting  at  Z  -  - 12  cm. 

The  temperature  gradient  at  Z  -  -14  cm,  corresponds  to  the 
spanwise  edge  of  the  film  cooling  for  blowing  rate  of  M  -  0.98. 

c.  Effect  of  spanwise  cortex  position. 

To  study  the  effect  of  the  vortex  position  relative  to  the  film 
cooling  injection,  the  vortex  position  was  changed  from  A  to  B  to  C  locations. 
Distributions  of  St/Sto  are  shewn  in  Figures  33,  37  and  40,  distributions  of 
St/Stf  in  Figures  34,  38.  41  and  distributions  of  T  -  Too  are  given  in  Figures 
35.  39  and  42.  Regardless  of  the  vortex  position,  the  following  phenomena 
were  observed  :(1)  on  the  upwash  side  of  the  vortex  a  low  heat  transfer 
region  exists,  where  St/Sto  is  less  than  for  a  film  cooled  turbulent  boundary 
layer  only,  (2)  on  the  downwash  side  of  the  vortex  a  wide  region  of  high 
St/Sto  exists,  because  of  thinning  erf  the  boundary  layer  by  the  vortex,  and 
(3)  at  X/d  -  7.4,  the  film  cooling  dominates  the  flow  and  the  vortex  seems  to 
have  very  little  effect  on  the  heat  transfer. 


Comparisons  of  Figures  33,  37  and  40,  show  that  tor  low  .wd 
values  up  to  ;7  3.  a  double  peak  in  the  St/Sto  distributions  is  present  for 
vortex  positions  A  and  B.  At  X/d  =  33.6  and  further  downstream  for  vortex 
positions  A  and  C,  the  double  peak  in  St/Sto  distributions  starts  to  grow 
spanwise  enlarging  the  portion  of  the  plate  covered  by  the  hot  spot  In  these 
cases  St/Sto  peaks  are  not  as  high  as  for  the  vortex  position  B  case  v-'gure 
331  A  maximum  St/Sto  of  1  i)3  was  observed  for  the  vortex  position  -■  .u 
X/d  =  %.b,  or  about  30  percent  greater  than  the  St/Sto  with  film  eooline 
only.  Low  St/Sto  and  St/Stf  regions  on  the  upwash  side  of  the  vortex  were 
observed  for  all  three  vortex  positions.  These  were  less  significant  for 
position  C,  while  varying  between  2-  6  percent  for  the  other  two  vortex 
positions. 

St/Stf  distributions  from  Figures  34,  33  and  40  show  that  the 
location,  shape  and  size  of  the  hot  spot  vary  depending  on  the  vortex  span- 
wise  position.  In  Figure  34,  for  vortex  position  B,  this  hot  spot  is  located 
from  the  middle  length  of  the  plate,  growing  in  size  as  the  vortex  develops 
downstream.  For  vortex  position  A,  the  hot  spot  is  shifted  upstream  and  is 
not  as  wide,  nor  as  long  as  for  position  B  Fot  vortex  position  C,  the  hot  spot 
begins  early  upstream  but  persists  further  downstream,  but  it  does  not 
reach  the  end  of  the  test  section,  and  it  is  very  narrow  A  region  of  high  heat 
transfer  is  observed  in  all  three  cases,  but  this  area  is  larger  for  vortex  posi¬ 
tion  A,  covering  about  2/3  of  the  test  plate.  A  heat  transfer  gradient  is 

1 

aligned  along  the  vortex  center  for  all  three  vortex  positions,  but  it  is 
steeper  and  more  persistent  for  vortex  position  B  The  region  of  low  heat 

I 

I 
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transfer,  where  Si/Stf  is  less  than  1,  is  largest  for  position  B  and  smallest  "or 
position  C 

Figures  35.  39  and  42  show  the  T  -  Too  temperature  field  in 
degrees  C  for  X/d  =  41.9.  For  these  results  the  heat  transfer  plate  was 
maintained  at  the  t'reestream  temnerature  and  the  coolant  was  heated  to  s  ; 
-C  i  unheated  plate  i  Overall  qualitative  trends  for  all  three  \urte\-  pus  it.  on/, 
are  similar  coolant  is  lilted  away  the  wail  bv  the  uowash  side  ot  the  ■  orev 
on  the  downwash  side  of  the  vortex,  secondary  flows  cause  the  freestreum 
fluid  to  be  located  very  near  the  wall,  making  the  boundary  layer  very  •.hir 
and  minimizing  the  effect  of  the  film  cooling.  A  significant  spanwise  gradient 
in  temperature  exists  which  is  located  near  the  vortex  core,  and  extends 
about  3  cm  from  the  wall.  As  expected,  the  spanwise  location  of  this  gradi¬ 
ent  changes  with  vortex  position.  The  effect  of  the  downwash  side  r,  the 
vor .ex  is  most  significant  for  vortex  position  B.  For  this  case  Figure  35a 
shows  that  the  boundary  layer  is  thinned  such  that  it  only  extends  2  cm 
from  the  wall.  Significant  changes  also  occur  near  the  upwash  side  of  the 
vortex  where  secondary  flows  from  the  vortex  convect  film  coolant  awav 
from  the  wall  The  most  significant  changes  in  this  part  of  the  fiow  occur  for 
vortex  position  C. 


Measurements  were  made  at  Uco  =15  m/s,  with  a  blowing  ratio 


M  of  0.47,  in  order  to  determine  the  effect  of  low  blowing  ratio  on  heat 
transfer.  These  tests  were  also  conducted  to  obtain  additional  verification  of 
general  conclusions  from  results  obtained  at  10  m/s  and  M  -  0.98  blowing 
ratio.  For  these  tests,  vortex  generator  n2,  was  used  at  positions  A,  B  and  C. 


Results  are  shown  in  Figures  nj  —  48.  Distributions  of  St/Sto  are  sinAvn  in 
Figures  4\  and  i7  and  distributions  of  St/S tf  are  given  in  Figures  ^  *<> 
and  48. 

From  a  general  qualitative  point  of  view,  the  results  at  15  m/s  and 
M  =  0  47  are  similar  to  those  of  10  m/s  and  M  =  0.98  The  most  sigrui  .cur : 
a  u  an  1 1  tut  me  uii. erence  is  mat  St/Sto  maxima  are  nigner  and  ^t/.'to  m  •- 
are  lower  u»r  Xi  =  0  h7  in  addition,  the  peaks  of  local  neat  transler  see m 
increase  more  rapidly  with  downstream  distance.  At  the  downstream  end  oi 
the  plate  where  X/d  =  96. o.  local  St/Sto  for  M  -  0.47  mav  be  as  high  a?  i  1  l 
compared  to  about  1.05  for  M  =  0.98.  The  rapid  downstream  growth  ot  these 
maxima  becomes  more  apparent  when  one  considers  the  film  cooled  bound¬ 
ary  laver  alone,  since  St/Sto  ratios  at  M  =  0.47  are  lower  than  those  at  M  = 
0  98. 

At  M  =  0.47,  the  effect  of  changes  in  the  spanwise  position  of  the 
vortex  are  much  less  significant  than  at  M  =  0.98.  At  higher  blowing  ratio 
many  quantitative  changes  in  the  local  heat  transfer  result  as  the  vortex 
position  is  changed.  In  particular,  primary  and  secondary  peaks  of  St/Sto 
were  observed  at  X/d  =  17.5  and  33  6  for  vortex  position  B.  and  X/d  =  54  o 
and  75  6  for  vortex  position  C.  At  the  M  =  0.47  blowing  ratio,  the  shapes  ot 
St/Sto  distributions  for  different  vortex  positions  are  very  similar,  especially 
for  X/d  --  54.6,  75  6  and  96.6.  Generally,  quantitative  changes  for  M  =  0  47 
occur  only  in  the  spanwise  locations  of  the  St/Sto  maxima  and  minima.  At 
X/d  =7  4,  the  effect  of  the  vortex  position  changes  is  minimal  and  the 
spanwise  heat  transfer  rates  are  mostly  affected  by  the  film  cooling  At  X/d 
-  I  7  5  both  primary  and  secondary  St/Sto  peaks  are  present  for  ail  three 


vortex  positions.  The  only  significant  St/Sto  changes  with  vertex  position  ,c- 


cur  at  X/d  -  33  6  where  a  small  secondary  St/S  to  peak  is  observed  when  the 

vortex  is  located  at  position  A  only.  Such  behavior  indicates  that  the  effects 
of  film  cooling  at  M  =  0.47  do  not  persist  as  far  downstream  as  at  higher 
blowing  ratios.  In  addition,  the  vortex  seems  to  more  completely  dominate 
the  flow  i  ield  as  the  blowing  ratio  decreases. 


4.  Free  stream  10  tn/s,  blowing  ratio.  M  -  1.26.  Fffect  of 

high  blowing  ratio. 

In  order  to  study  the  effect  Oi  n«gh  blowing  ratio  on  the  film  cooled 
turbulent  boundary  layer  with  embedded  vortex,  measurements  were  made 
at  Uoo  =  10  m/s.  with  M  =  1.26  and  the  vortex  at  position  B.  Results  are 
shown  in  Figures  49-  52. 

The  St/Sto  and  St/Stf  data  in  49  and  50  show  the  same  overall 
qualitative  trends  shown  by  measurements  at  M  *  0.98,  0.47.  and  0.6 8 
augmented  heat  transfer  near  the  downwash  side  of  the  vortex  and  reduced 
heat  transfer  near  the  upwash  side.  Also,  at  X/d  »  7.4,  film  cooling  rather 
than  the  vortex  dominates  the  spanwise  variation  of  heat  transfer.  The  most 
important  quantitative  differences  between  there  results  and  those  at  other 
blowing  ratios  are  the  magnitudes  of  St/Sto  with  film  cooling  only.  Figure  26 
shows  M  =  1.26  St/Sto  data  at  10  m/s  which  is,  on  the  whole,  greater  than 
data  at  M  =  0.68,  but  just  lower  than  results  for  M  =  0.98.  Results  presented 
in  Figure  49  reflect  the  same  overall  trends. 

Details  of  the  spatial  variations  of  normalized  Stanton  numbers  in 
Figures  49  and  50  are  now  discussed.  At  X/d  =  17.5,  33  6  and  54.6  double 
peaks  in  St/Sto  distributions  are  present.  These  seem  to  result  from  inter- 


action  between  cooling  jets  and  the  vortex,  occurring  at  Z  -  -1 — 3  cm  ana  at 

7.  =  2  -  3  cm. 

Contours  of  T  -  Too  l^C'l,  at  X/d  =  41.9  in  Figure  Si  show  that  the 
peak  near  *2-3  cm  results  from  the  downwash  of  tree  stream  fluid  to  near 
wall  regions  bv  the  vortex.  The  corresponding  St/Sto  peak  increases  in 
magnitude  as  the  the  boundary  iaver  cunvects  downstream  reaching  mag¬ 
nitudes  as  high  as  i .0 5  at  X/d  75.6  and  96.o.  At  these  downstream  locations 
augmented  St/Sto  cover  a  large  spanwise  portion  of  the  film  cooled  test  sur¬ 
face  extending  from  -5  cm  to  O  cm  (Figures  49  and  5df  Magnitudes  ot 
St/Sto  maxima  in  these  hot  spots  at  M  =  126  are  about  the  same  as  those  at 
M  =  0.98. 

Figure  50  shows  the  heat  transfer  rate  over  the  test  plate  in  terms 
of  St/Stf.  Here,  the  hot  spot  is  located  from  X/d  =  54.6  to  X/d  =  96  6.  This 
region  is  very  large,  confined  to  the  downwash  side  of  the  vortex.  The  tem¬ 
perature  gradient  on  the  upwash  side  of  the  vortex  is  not  very  steep  in  some 
streamwise  locations.  A  comparison  between  Figure  50  and  Figure  32  where 
M  =  0.68,  show  that,  the  hot  spot  is  present  at  the  same  spanwise  and 
streamwise  location  in  both  cases,  but  in  50  the  hot  spot  is  thinnei ,  and  the 
area  of  low  St/Sto  on  the  upwash  side  of  the  vortex  is  larger 

Temperature  contours  at  X/d  =  419  in  Figure  51,  show  that  the 
St/Sto  peak  at  7  =  -2—3  cm.  lies  beneath  a  region  of  large  temperature  gra¬ 
dients  very  near  the  vortex  center.  Figure  52d.  shows  strong  secondary- 
nows  at  this  location  which  apparently  sweep  coolant  further  from  the  vor¬ 
tex  center  f to  smaller  Z),  such  that  it  collects  near  -5  Z  -8  cm  Contours  of 
vorticitv  magnitude,  calculated  from  these  secondary  flow  vectors  for  \1  - 


1.26  vortex  position  B,  are  shown  in  Figure  52e.  These  show  that  (1)  the 
core  of  the  primary  vortex  is  located  at  Z  -  -2  cm.,  Y  =  3.5  cm.,  and  (2)  that 
the  secondary  peak  in  heat  transfer  is  not  a  result  of  any  form  of  secondary 
vortex.  Augmented  mixing  is  initially  created  by  the  secondary  flows 
beneath  the  vortex  center,  resulting  in  locally  higher  heat  transfer.  This  se¬ 
condary  peak  shown  in  Figure  49  at  X/d  =  17.5,  33.6  and  54.6  eventually 
becomes  indistinguishable  at  larger  X/d  equal  to  75.6  and  96.6.  The  se¬ 
condary  peak  is  bigger  and  more  persistent  than  at  M  =  0.98,  indicating  that 
the  shear  layer  interaction  causing  higher  heat  transfer  is  dependent  upon 
the  mass  flux  of  film  cooling. 

Another  important  feature  of  the  results  in  Figure  5 1  is  the  gradient 
of  heal  transfer  present  at  Z  -  -10-  -12  cm.  for  0  <  Y  <  5.0  cm.  This  gradient 
is  present  at  the  edge  of  the  film  cooled  region  :  coolant  was  injected  from 
only  7  holes  extending  across  a  span  from  -9  cm.  to  +9  cm. 

Additional  flow  field  information  at  X  -  1.48  m,  X/d  -  41.9  is  given 
in  Figures  52a.-c  for  10  m/s,  M  -  1.26  and  vortex  position  B.  Data  was  taken 
using  a  five-hole  pressure  probe,  [Ref.  1.:  pp.  14-151,  to  measure  pressure  at 
800  points  in  a  spanwise  plane  in  order  to  further  investigate  the  effect  of 
the  high  blowing  ratio  in  the  turbulent  boundary  layer  with  embedded  vor- 
tei.  These  results  are  consistent  with  those  of  Evans,  [Ref.  1 1.  They  show  :  ( 1 ) 
that  high  velocity,  high  total  pressure  fluid  is  swept  near  the  wall  on  the 
downwash  side  of  the  vortex,  (2)  a  velocity  and  total  pressure  deficits  exist 
at  the  center  of  the  vortex,  and  (3)  low  velocity  fluid  is  swept  away  from  the 
wall  near  the  upwash  side  of  the  vortex.  Vorticity  contours  in  52e.  show  a 


large  region  of  negative  vorticity  in  the  upwash  region, (-8  <  Z  <  -4  cm).  In 
addition,  two  other  regions  of  negative  vorticity  result  from  the  secondary 
flows  present. 
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IV.  SUMMARY  AND  CONCLUSIONS. 


Heat  transfer  measurements  were  made  of  turbulent  boundary  layers 
with  film  cooling  and  embedded  vortices  for  freestream  velocities  of  10  and 
!  15  m/s  .  To  obtain  the  data,  a  heat  transfer  test  surface  was  designed  and 

developed  to  provide  constant  heat  flux  over  its  area,  with  126  embedded 
thermocouples  for  detailed  spanwise  resolution  of  temperature.  Extensive 
qualification  tests  show  that  the  surface  gives  excellent  spatially  resolved 
heat  transfer  coefficients  over  an  area  of  43.8 1 5  cm  x  1 1 1 1.76  cm. 

Baseline  measurements  show  excellent  agreement  with  Stanton  number 
correlations  for  a  flat  plate  with  constant  wall  heat  flui  and  unheated  start¬ 
ing  length.  Results  of  turbulent  boundary  layer  with  embedded  vortex  show 
excellent  agreement  with  Joseph,  (Ref.  3-1.  and  with  the  literature.  Results  of 
turbulent  boundary  layer  with  film  cooling  at  different  blowing  ratios  show 
|  expected  trends,  consistent  with  the  data  of  Goldstein  and  Yoshida,  [Ref.  18) 

Longitudinal  vortices  cause  significant  changes  in  heat  transfer  and 
structural  characteristics  of  film  cooled  turbulent  boundary  layers  Heat 
transfer  augmentations  as  large  as  30  percent,  and  reductions  as  high  as  10 
t  percent  were  observed  to  persist  as  many  as  23  boundary  layer  thicknesses 

or  96  film  cooling  hole  diameters  downstream  of  film  cooling  injection  loca¬ 
tions.  The  effects  of  the  embedded  vortex  on  heat  transfer  in  film  cooled 
boundary  layers  are  significant  and  important : 

1.)  Near  the  downwash  side  of  the  vortex  the  heat  transfer  is  augmented, 
vortex  effects  totally  dominate  the  flow  behavior,  and  the  effects  of 


$ 
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the  film  cooling  are  negated.  '*  Hot  spots "  will  exist  for  blowing  ratios 
ranging  from  0.47  to  1.26,  and  vortex  circulation  to  freestream  velocity 
ratios  of  about  1.6  cm. 

2. )  Near  the  upwash  side  of  the  vortex,  the  coolant  is  lifted  off  the  wail  and 

pushed  to  the  side  of  the  vortex,  increasing  the  surface  area  protected 
by  the  film  cooling. 

3. )  Changing  the  position  of  the  vortex  with  respect  to  the  film  cooling  jets 

results  in  significant  local  quantitative  changes  in  heat  transfer  occur, 
even  though  the  overall  qualitative  trends  remain  unchanged. 

4. )  Near  film  cooling  injections  locations  for  X/d  up  to  7.4,  the  film  cooling 

dominates  the  flow  behavior  and  the  vortex  seems  to  have  very  little 
effect  on  spanwise  variations  of  heat  transfer. 

Results  (1)  and  (2)  are  mostly  a  consequence  of  the  intense  secondary 
flows  produced  in  the  plane  perpendicular  to  the  axis  of  mean  vorticity. 
These  results  are  consistent  with  those  obtained  by  Joseph,  (Ref.  3.:  pp.  54], 
and  Evans,  (Ref.  1.:  pp.  28 j. 

At  high  blowing  ratios: 

1. )  A  double  peak  in  the  St/Sto  distributions  was  observed  to  occur 

between  X/d  -  7.4  and  54.6. 

2. )  The  change  in  spanwise  position  of  the  vortex  in  relation  with  the  film 

cooling  jets  affects  the  magnitude,  shape  and  spanwise  position  of 
St/Sto  peaks,  and 

3  )  Secondary  St/Sto  peaks  become  higher  in  magnitude  and  more 
persistent  with  respect  to  downstream  development  as  the  blowing 
ratio  increases.  The  double  peaks  observed  at  M  -  1.26  were  not  due 
to  a  secondary  vortex,  but  to  an  interaction  between  the  vortex  and 
film  cooling  which  depends  on  blowing  ratio. 


At  low  blowing  ratios: 

1. )  The  change  in  spanwise  position  of  the  vortex  has  very  little  effect  on 

local  heat  transfer  distributions  except  to  change  the  locations  of 

Stanton  number  minima  and  maxima,  and 

2. )  St/Sto  distributions  exhibit  only  one  peak  which  increases  in  magni¬ 

tude  with  downstream  development. 

It  is  recommended  that  flow  visualization  study  of  the  interaction  of  the 
vortex  and  film  cooling  be  conducted  in  order  to  enhance  the  understanding 
of  some  of  the  complex  phenomena  observed  during  the  course  of  this 
study. 


Fig.  2.  Vortex  Generator  Position. 
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Fig.  33b.  Spanwise  Variations  of  Stanton  Number  Ratios. 
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FREE  STREAM  10  M/S,  WITH  FILM  COOLING,  VORTEX  GEN. *2  AT  4 . 79  cm 

Fig.  33e.  Spanwise  Variations  of  Stanton  Number  Ratios. 
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Fig.  40b.  Spanwise  Variations  of  Stanton  Number  Ratios. 
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Fig.  40d.  Spanwise  Variations  of  Stanton  Number  Ratios. 
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FREE  STREAM  15  M/S,  WITH  FILM  COOLING,  VORTEX  GEN.  AT 

Fig.  45c.  Sparwise  Variations  of  Stanton  Number  Ratios. 
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APPENDIX  3 


UNCERTAINTY  ANALYSIS 

L  near tain tv  analysis  was  periurmed  using  uie  met aou  or  wc-si;.;  .iv 
kline  and  McGmtnck  [Ref.  19|.  This  is  the  rooi  sum  square  me  thou,  where 
the  uncertainty.  <5p  .  of  some  function  F,  is  a  function  of  the  independent 
variables  Xn,  according  to  the  following  expression; 


All  uncertainties  are  based  on  9 5  confidence  levels.  To  calculate  the 
uncertainty  of  the  Stanton  number,  it  is  necessary  to  calculate  first  the 
uncertainty  in  the  heat  transfer  coefficient,  h.  The  following  independent 
\  ariable  uncertainties  were  determined  :  5q  =  t  16  \V/m2.  5tw  =  -  0.5  20.  $t.v 
-  ?  0  ;  -C.  Here,  the  uncertainty  of  heat  ioss  by  convection  was  based  on  a 
3%  error  in  radiation  losses  and  2.5%  error  in  conduction  losses.  The 
uncertainty  of  Tw  is  higher  than  Too  due  to  higher  uncertainty  in  the 
calculation  of  the  contact  resistance.  From  these  parameters  the  uncertainty 
of  h  was  determined  to  be  4.8%  or  approximately  ±  i  7  W/m<-  WC,  based  on  an 
h  value  of  35  W/m2^. 

To  calculate  the  uncertaintv  of  the  Stanton  number  besides  the  uncer¬ 
tainty  of  the  heat  transfer  coefficient,  the  following  independent  variable 
uncertainties  were  determine  ;  bpoo  =  r  0^1  :■  *:j/ rn-%,  Suoo  =  zQ.5  uVs  = 
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=  •  -  -■  -  -  -'K.  The  uncertainly  of  Cp  was  based  on  the  assumption  <.a  constant 
properties.  From  these  parameters  the  uncertainty  of  the  Sian :•••;*.  ..  me; 

was  calculated  to  be  t  6.04;",  or  approximately  1.87  x  i(M  based  on  a  typical 
Stanton  number  of  0.0035-  The  uncertainty  of  St/Sto  is  derived  from  the 
parameters  already  mentioned  and  was  estimated  to  be  :  ~  .17  or  :  o'v 
based  tm  a  St/Sco  of  0.75. 

To  calculate  the  uncertainty  of  the  blowing  ratio.  VI.  along  die  uncer¬ 
tainty  values  determined  for  Uoo  and  poo.  the  following  independent 
variable  uncertainties  were:  5i;jnj  -  t  i.C  m/s.  and  5pjnj  =  ;  cl  :  "he 

uncertainty  of  U1Qj  is  higher  than  Uoo  due  to  higher  uncertainty  of  the 
rotometer  and  the  area  of  the  injection  holes,  and  the  uncertainty  of  pjnj  is 
also  higher  than  pco  due  to  higher  uncertainty  in  the  estimation  of  Tmj  from 
Tp.  From  these  parameters  the  uncertainty  of  M  was  calculated  to  be  -  8 
or  about  t  0  TCTo  for  a  M  value  of  0.98. 


j 

< 
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APPENDIX  C 


SOFTWARE 


ne  iouv Avmg  programs  are  .isied 

eOci'gv  OUiuuCe  estimation  j  of  COituUdioil  iOS.'c: 


; 1  •  ;  r,  r. 
« .  •  l\  A  *  > 


y  vt :r  i 

.1.  i L  ‘  ,  l 


heal  transfer  program  lo  acquire  thermocouple 
readings  from  D.\S..  and  store  information  in  data 
files  in  floppy  disk 


STANFC2  :  heal  transfer  program  to  calculate  Stanton  numbers 
for  all  flow  conditions 


AG/fPRO  :  neat  transfer  program  to  acquire  temperatures  from 
automatic  traversing  device  for  temperature  surveys 

PLOTRUN  :  program  to  read  a  data  file  to  plot  temperature 
contours 

PTSTAV  :  programs  read  a  data  file  to  plot  spanwise  averaged 
Stanton  numbers  vs.  Reynolds  number 


PTSTLC  :  reads  a  data  file  to  plot  spanwise  local  heat  transfer 
coefficients 


PLSTRTIO  :  plots  spanwise  St/Sto  for  film  cooling  only. 

PLSTRVOR:  piots  spanwise  variations  of  St/Sto  for  embedded 
vortex  data  onlv 


PLS7RYV.  plots  spanwise  variations  of  St/Sto  for  film  coding 
only  and  St/Sto  for  film  cooling  and  embedded  vort 
bv  rows. 
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5URFC0NT:  plots  surface  contours  of  St/Sif 


PL ST RFC 


plots  spanwise  variations  of  local  St/Sto  for  i 
cooling  data  only 


Name  of  variables  used  are  intended  to  be  self  explanatory 
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i  PROGRAM  ENERGE 

i 

! THIS  PROGRAM  ACQUIRES  MULTIPLE  CHANEL  THERMOCOUPLE  CATA 
! AND  PERFORMS  ENERGY  BALANCE  TO  ESTIMATE  CONDUCTION  LOSSES 
ILI5RANI /ORTIZ  VERSION,  JUNE  1 ScB 

i 

DIM  £( 238  ;  ,T( 233  ) 

! CHANNELS  3-73  ,  i 33- i 45 
! COPPER  CONSTANT AN  THERMOCOUPLES 


1  2  2 

rnfu^M"  *  •*  » 

rRi.ncr  -  3  .  tj  • 

•  jw 

P  R  i  N  t 

) 

PRINT  " ENERGY  BA 

152 

PRINT  + 

1 S0 

!  ENTER  POWER  IN 

:  S  : 

DISP  “  i HIT , CONTI 

1  5  2 

-A  f\  .  i  r\  • 

rnu’cc 

\  72 

DISP  "  cNTER  SUP  Pi 

130 

INPUT  ro p  i  Uc  Its 

1  30 

PRINT 

200 

PRINT  "CURRENT  = 

2  » G 

PRINT 

•—  'ii 

FRINT  "TEMPERATU! 

23  3 

PRINT 

2+3 

Td'vS-3  . 

Z50 

PRINT 

2S3 

.“OR  I  -  i  TO  i  3 

273 

OUTPUT  735 ; " A I " s 

233 

ENTER  703 iK 

2  gw 

EC  I  /=.'<*;  333332  , 

330 

T(  I  >.01 SZO5+.0Z: 

3  5  C 

Tave-Tave+T ' I  ! 

3Z3 

PRINT  USING  333; 

333 

IMAGE  DOC  ,-■/  ,DDC1 

34-3 

NEXT  I 

3  S3 

PRINT  "+*♦»»**** 

36  3 

PRINT 

«J  i  w 

PRINT 

333 

FOR  1-43  TO  SS 

3S3 

OUTPUT  705; “ ft I “ ; 

*  r*  r* 

ENTER  7  3  3  ;  X 

4  i  r» 

4.  1  E 

E‘!  I  /-X  *  i  3203G3  . 

42G 

T <  I  >. 013235a.  32! 

430 

Tave“Tave+T( I > 

440 

PRINT  USING  453; 

453 

IMAGE  ODD, ZX, DCS! 

462 

NEXT  I 

473 

PRINT  '•*-**-»*♦* 

4  8  G 

PRINT 

4S0 

Ta  va-T avs.‘4S . 

523 

PRINT  "AVERAGE  FI 

LY  CURRENT  IN  AMP'S,  ANC  VOLTAGE  IN  VCl T~‘ 


MEASURED  ( CE5  C> 


;  T  a  v  e 


.r.v.v.v.v.v  v.v  'wiK  rjrrji 


540 

PRINT 

550 

'GET  T  AMBIENT 

560 

OUTPUT  73S;"AI"i3;" 

VT !  '• 

570 

ENTER  709 s X 

530 

£(  3 )*X* 1 330 

530 

T( 3  1=25  .573 ♦£(  0  )  - 1  . 

936879*E(3'-e/sl+ 

)*Ev 

•3  )*E(  C  i 

6  1 0 

Tarcb^T ( 0  ) 

S20 

PRINT  " AMS  I ENT  TEMP 

.  = " ; Tarcb 

630 

PRINT 

S40 

Tax ff=Ta^e-Tanc 

852 

PRINT  "TPLATE-TAM9 

= "  i T  cU  f  f 

652 

F  F I N  T 

£  -**  /\ 

•4  <  W 

' SET  AV ERASE  TEMP  I 

N  I wSbLA i I  ON 

bow 

T l ns  1 =3 . 

53G 

FOR  1=147  TG  143 

703 

OUTPUT  73S > ’AI " i I ; " 

VT  1  " 

710 

ENTER  739 iX 

T"* 

E  (  I )  =  X  *  1 332  . 

.-1  •—  r\ 

T  (  I  > = 2  5 . 5  T  3  »  E  \  I  1-1  . 

32-6S79*E(  I  ')*•£(  I  )- 

)  *E \ 

I  >*E<  I  / 

750 

Tin5iaTinsl+T( I > 

763 

'PRINT  USING  4c3 ; I  , 

El  I  ) ,T( I  ) 

«-•  •-»  r% 

l  l  W 

NEXT  I 

733 

Tins! aT ins./ 3. 

790 

Tins2*0 . 

£33 

FOR  1=153  TO  152 

3 !  0 

OUTPUT  733 i "AI" s I ; " 

VT  1  " 

i- .251 


’  *  £  (  3  1  *■  E  {  3 


323  ENTER  705 ;X 
332  Ev  I  !«X*  \  000  . 

343  7  ( I  1=26 . 573  *£■:  I  )-|  .33S8?3*E>  I  )*E(  I  >-,9378S*E(  I  >*E(  I  >*£■:  I  .  26  1  Z77*E<  I  )*£i  1 

>*£<  I  :■*£<  I  > 

850  'PRINT  USING  450; I  ,E( I > ,T< I ) 

370  Tins2*Tins2  +  T( I  / 

383  NEXT  I 

8SS  Tir,52  =  Tins2/'3. 

300  PRINT 

3 1 0  PRINT  "TINS  1  =";Tins1  TIMS2  -"jTins2 
3  2  G  PRINT 

933  'ENTER  VALUES  GF  THERMALC0N0UC7. ,  AREA,  AND  THICKNESS  OF  INSULATION 
943  K  => .  3 4  1  W/M  0E3  C 

353  Da -.32  54  !  M 

563  A=.4337  ' M " 2 

3<0  iOALCULATt  HEAT  FLUX  THROUGH  INSULATION 
3S3  Gins=K»A*ITinsi-Tins2 )/Da 
330  'GET  Q  IN  ( WATTS  ) 

1000  Qi n=Vol t s ♦Anps 
1213  Qcond=Q;n-Qlns 

1023  PRINT  "G  IN  =  ";Gir. GIMSULATIGN  •NQins,"  GGONDUCT ION  =“:Gcona 
1033  PRINT  ’♦•***»»****»****»**»*»***»****«»4*»»»«******»».*'- 

1340  STOP 
! 052  END 


169 


in  m  m  m 


!  nt.i.'IOLiOUPLc  REnuiNG- 


!  RLguiRC  :  ncn.'lUL UurLc 
PRINT 

FOR  I *  1  TO  73 
OUTPUT  7  3  3 ;  "  A I  “  1 1 ;  “  V  7  i  " 
ENTER  7G3;X 
E  i  I  =X  •  i  330*380 


BOO 

T(  I ' “ . 3 ! ,3205+ . 3Z5G4S »E ’  I 

/- . 0202C3S3 1 *E(  I  1*E\  I  ) 

s  ;o 

PRINT  USING  523 i I  ,T( I  } 

520 

IMAGE  GDDt4X,0Q.CD 

500 

OUTPUT  SFaih_!|7iI) 

5  AO 

NEXT  I 

353 

1 

3  50 

-CP  I » i £  3  T C  ‘ A  5 

<-»  •**  r\ 

C  *-J 

OUTPUT  70S : “ A I " ; I ; " UT ’  " 

UW  V 

ENTER  703  *  X 

£33 

E-:  I  ;*:<-■  300030 

700 

Ti  I i-.3t 3Z5-.22S54S-E-  I  ; 

-.030000531  +E-:  I  :•*£•:  I  > 

7 » c 

r  -  n 

■'  <J  'c' 

OUTPUT  3Fath_i j  T' I  ) 

"A3 

NEXT  I 

'GET  AMBIENT  TEMPERATURE 

?S£- 

•0 UTPUT  7 3 3 : "  A I "  j S  ;  “ 'J T  • " 

<■  i  £ 

ENTER  723  ;X 

753 

E< 2  /“X* ! 030 . 

"1  ,n  n 

{  ^  0 

T i  3  ) ■  2 5  .  S 7 ~ E  i  2  > - 1  .33537 

5  ♦  E  \  3  :  +  z.\  0  f r .  3  3  7  9  5  ♦  £  i  3  }  * 

c  *2  3 

T  chanfc*  T x 0 > 

313 

PRINT  “ THERMOCOUPLE  A ME  I 

ENT  ! EM? . ( DEG  Cl  =  ‘ «  Tcha 

620 

PRINT 

830 

'GET  THE  REST  OF  THESHOC 

CUPLES 

543 

FOR  I*; 47  TO  ISO 

350 

OUTPUT  70S i " AI " ; I s "UT1 " 

©SO 

ENTER  73S»X 

373 

£■■  I  1-X+13CS. 

/n  1 

o  ■:«; 

T ;  I  ’ZB  .  573  »E I  i  - 1  .  33537 

3*E<  I  )*E'  I  /  + .  327375-E(  I  1 

530 

f3  R  I N  T  USING  0  2  3  ;  I  ,  T  ■'  I  ) 

500 

OUTPUT  OF  at  h_!  :  T  •!  I  } 

3  i  0 

ME  X  7  I 

303 

ASSIGN  <?Fath_l  TO  * 

300 

i TRANSFORM  TO  SI  UNITS 

340 

i 

350 

R.amb3? an-.b*33c5 . 22 

1  PASCALS N  i  M  '  2  ) 

363 

Oeltao=Dsi tap*24S.7 

! Pascals;  N/r-r  z  :• 

373 

Fstanp^Ti i 47  >  +  273. IE 

*DtG  KELU I N 

330 

Eo“Parb / ( 237  *Fs  t  ercp  ) 

•’AIR  DENSITY  (KG/M3! 

930 

U  i  n f  » <  2  » Oe  1 1  ap  / Ro  i  “  ,  5 

1  FREE  STREAM  VELOCITY 

1  G30 

[ 

io;o 

PRINT 

I  003 

OUTPUT  OPa  t  h_2  ;  Ru.nnc  ,  Anp 

-j  , Volta  , Parch  , Oe  i  t  ap  ,  Tarcb 

1030 

ASSIGN  3Pa th_2  TO  ♦ 

i  '*•  \J 

l 

!  253 

GOTO  ! 1  S3 

issa 

’SUBROUTINE  TO  CALCULATE 

FILM  C OGLING  Dnin 

:-E(o;-i 


Ji nr  ,Ro  , Fs tercp , An 


.'TOWWVITOTO ^jv.-'#.  -v#w  -wu<  rj*rj*  r*  -ji  -  ,• 


jj 

I 


t 


f 


i 

I 

N 


1070  PRINTER  IS  1 

I0S0  PRINT  "ENTER  PLOW  %  FROM  RGT0ME7ER " 

1030  INPUT  r  I  w 

}\M  PRINT  "ENTER  PLENUM  FRES3URE  DIFFERENCE  i IN  H20 ) “ 

I  i,l  0  INPUT  Oelpi 
1123  Qelpi”0eipi*248.7 
1 i 33  Flw=Fiu* .303033455 
1140  CREATE  SCAT  "FILQT2S ' ,2 ,3 

HS3  ASSIGN  3Path3  TO  "FILDT2qH  ' OPEN  FILE  TC  RETRIEVE  FILMCCCLING  DATA 
1163  OUTPUT  $Paih3;Fiw,D©Jpi 
1173  ASSIGN  @Path3  TO  * 
liae  GOTO  1 93 

1133  MASS  STORAGE  IS  INTERNAL  ,4" 

i  Z33  end 


p 

p 


I 

s 

\ 

s 

\ 


> 
t 
t t 

i 


> 

1 


172 


f  VVV.VVW^V.«V.>\*V/  UV/V' 


>  .\  .V, 


.V 


's’J’-Jtrx'J'-J*  v  v -/W r> -> -> ->  V V  V  ->  ■>■ 


r 


i  k 


10 

23 

30 

40 

50 

S3 

7G 

80 

30 


FRCSRAM  STANFCI 


!  1 0  .1/3  FREE  STREAM  VELOCITY,  13  INJECTION  HOLES  OFEN 
?  THIS  PROGRAM  ACQUIRES  MULTIPLE  CHANEL  THERMOCOUPLE  DATA 
1  CALCULATES  HEAT  TRANSFER  COEFFICIENTS  AND  STANTON  NUMSEPS  RATIOS 
! AND  FILMCGGLINS  PARAMETERS. 

! LIGRANI/ORTI2/ JOSEPH  VERSION,  NOVEMBER  1335 
i UPDATED  BY  ORTIZ,  JULY  1987. 


110  MASS  STORAGE  IS  *: INTERNAL  ,4 , 1' 
123  CREATE  BDAT  "STRFCV10" ,504,8 

120  ASSIGN  OPath  4  TO  “STRFCV10" 


*  *12' 

DIM  T 

(  230  •’  , M •  230 

1  ,  S 1 1  2  3  0  1  ,  X  !  (  5  ) ,  H 1  \  a  / ,  S 1 1  \  3  )  ,  X  2 1 

1  c 

*  20  4.' 

DIM  H 

o  \  2  3  3  ) ,  3 1 0  v 

203)  , X 1 0 ( 3 )  ,X2a<  200 ) :Zo<220) ,3tr 

153 

1 

1  70 

PRINT 

crir  2  d  ] 

4  /-«  O. 

1  OW 

t 

1  JiL’ 

r  R I  '1  i 

"rRE  YOU  {J 

SING  VORTEX  GENERATOR?  C  3 : Y  OR  4 

*■'  2  2 

i  ;N*r  “J  t 

V 

?  !  0 

IF  '{ = 

3  THEN 

21:0 

r*  'r. 

INT  H ENTER 

TyPE  OF  VORTEX  GENERATOR" 

230 

INPUT  Vert 

240 

PR 

INT  "ENTER 

LOCATION  FROM  CENTER  LINE  IN  cm" 

2  w£ 

IN 

PUT  L  0  c 

npii 

4.5  v 

ELSE 

->  —r\ 

A.  •  W 

■SOTO 

283 

2  S3 

END  I 

17 

230 

PR  INT 

ER  IS  731 

303 

PRINT 

CHRS <  2? )  i CHR$(  33  > ?  CHR3 !  03  >  s  CHR*<  S3  i  5 CHRS ( 

313 

PRINT 

CHRS ! 27  )  1  CHR3( 38  )  ; CHR  S ( I 35  )  : CHRS ( 43 ) ; CHRS ( 

OPEN  FILE  TO  RETRIEVE  STANTON  Nos  RATH 

;;  .z-zos  '  .St  r  230  . 


323 

f  f% 

343 

353 

353 

370 


1 


ASSIGN  OPa  th 


!  I  u  I 


OAFC 1 0" 


PGR  1=1  TO  133 
ENTER  aPathliTd) 

NEXT  I 

ASSIGN  SPathZ  TO  "IDAFCiO" 


s 

33G 

ENTER  3Path2;R 

't 

393 

PRINT 

» 

430 

; 

s 

11 

412 

DISP  "(HIT 'CON 

1 

423 

PAUSE 

% 

433 

IP  Arts*  i  THEN. 

- 

4. 

443 

PRINT  "STANTON 

- 

453 

ELSE 

453 

PRINT  "STANTON 

470 

END  IF 

1 

430 

PRINT  “******* 

1 

430 

PRINT 

/ 

500 

PRINT  "RUN  $  ( 

, 

510 

PRINT 

'' 

520 

IF  Y=3  THEN 

.temp  ,  Arts 


1***4H**»**t* 


173 


■  "V  /  V  "fc*  ■  „* 

v  \vV  O 


•% 


ra 

r, 


V 
v" 

V 

V 

vl 


i! 


A 


* 


A 


A 


3 


s 

v\ 

•a 


31 

$ 

.  j 


V 


530 

PRINT 

•VORTEX 

543 

PRINT 

552 

cLSc 

SS3 

GOTO  ESQ 

570 

ENC  IF 

530 

IF  Ana 

= !  THEN 

530 

1  PRINT 

OUT  DAT 

533 

PRINT 

'•AMPS  =“ 

S 1  3 

PRINT 

SZ0 

PRINT 

"AMBIENT 

533 

PRINT 

643 

1 

553 

PRINT 

“AIR  DEN 

553 

PRINT 

570 

PRINT 

"FREE  5T 

580 

PRINT 

533 

PRINT 

.*  rl  c  n .  :  _* .  ; : 

730 

3  PINT 

7;-3 

!  C.ALOU 

LATE  THE 

733 

■  s  e  -  C 

, 

i  **  u 

r*  r\r*  » 

P  J  M.  *  - 

i  •  W  •  w  w 

1  -a  #  ^ 

*  1  dvi*!  • 

753 

NEXT  I 

777 

7  a v  e=T 

dv6  !  4_ w 

print 

”  Av'EF^GE 

/  2  %j 

PR  I  NT 

330 

lENEr 3 

f  cnunNL; 

313 

t 

7  d  i  f  f  = 

■  av4* 1  •  j  h < 

3  J.C 

3cond* 

.  o53t  . 

n  . 

O-t/' 

1 dvota 

*  T  a  v  e  +  £  7 ! 

C2\j 

Tchafca 

=  Tchano<-; 

r*»~  -» 

ocw 

Grad* . : 

3C0G2S03 

"E  TEMPERATURE ,  MEASURED  (DEG  Cl  *"  (Tave 


Gi  n*nnpa*Vo i t 5 
Gcpnv*Q  >  n-Qcond-’Grad 
Oe 1 1  =7a  vaa  3  -Fs  t  ernp 


3  r  e 

PRINT  “7  PLATE  -  7 

333 

PRINT 

930 

PRINT  “ CONDUCTION  L 

543 

PRINT 

350 

PRINT  "RAD I AT I ON  LO 

950 

PRINT 

970 

i WALL  TEMPERATURES 

553 

1 

533 

Cp  * 1 005 .  1 

i  000 

Cr» .014  ! 

1010 

Ot*Cr*Qconv 

i  333 

G=Qconv/ . 4397  1 

1 033 

1 

1  040 

Tcavbs* 1 flvabj'Dt 

1353 

Oe i  t c  =  Tca»b5 -Fa  t  er-'p 

1350 

T he  t  a=  (  T 1  >  3b  -i -F 5 1  an 

IN  ( «iA TT  3  ) a '  i  G  i  n 


'  CEL7A  TEMP.  DUE-TQ  CONTACT  RESISTANCE 


1510 

Z  (  M  )  *  u  1  i  1 

t  520 

X2( M  )  =  X3 

1533 

FOR  I  =«-*-;  TG  J  *2  ! 

1S40 

X2( I >-X3 

1550 

i  >=z;  i-i  )+i  .27 

1550 

NEXT  I 

( 573 

X3=X3f . 2 

1630 

NEXT  J 

1630 

FOR  1  =  1  TG  125 

1730 

OUTPUT  @Fatn_4;Str< I ) ,Stfr< I 1 ,X2<  I) ,Z{  I 5 

1713 

PRINT  USING  1 720 1 1  ,  X2( I >  ,Z( I ) ,H( I  )  ,Str< I >  ,Stfr( I  ) 

1  7Z2 

IMAGE  DDD  ,2X ,DC.CD ,3X ,3D.DD ,5X .SQ.QQQE  ,5X  .SO.QDGE 

,5a  ,D . 3DE 

1  730 

NEXT  I 

1  7  43 

ASSIGN  i2r-ath_4  TO  *  i  CLOSE  Fli-r  WITH  STANTON 

NUMBERS  RATIOS 

1 753 

• 

1 CnLCULAT-  LOCAL  REYNOLDS  NUMBER 

i  ( V  u 

1 

'  733 

PRINT 

1753 

i 

Mu*. 33001  S3  ! CINEMATIC  VISCOSITY  FOR 

AIR  (M2/S1 

1313 

Faj*Uinf/ Nu 

13-0 

x  ■» ;  i  = 

1  833 

X  1 ' 2  '  = 1 .25 

1  34.0 

X 1 ( 3  )  -  1  .43 

1  353 

X 1 \ 4  > -  1 . S3 

1353 

X  1  1  5  )  =  1  . 3  0 

1373 

X 11  5  >  =  2 . 3  3 

1353 

! 

1333 

FOR  1=1  TO  5 

1303 

Rsv  v  I  :  =  1  ;  I  /*Fac 

1310 

NEXT  I 

l  -3-  u 

i 

1  333 

i CALCULATE  AVERAGE  STANTON  NUMBER 

1340 

FOR  I=i  TO  S 

1350 

H:  !• )  =0 . 

1953 

M=I *2 1 “23 

1  370 

FOR  J =M  TQ  1+21 

1  330 

HI  i  I  }=H  1  I  )tH(  J  ) 

1  533 

NEXT  J 

2033 

H 1  (  I  >=H 1 ( I  )/Z 1 

2310 

Still  >  — H 1  <  I  .'/ <■ Ro+Uinf  *Cp  ) 

2023 

NEXT  I 

2030 

PRINT 

2043 

i 

2053 

ASSIGN  @Fath7  TQ  “ST AU2 “ 

■*  *’*  c*  r* 
(.■Ob  0 

FOR  1=1  TO  5 

2370 

ENTER  3Ra t T,7 ; X 1  o(  I  )  ,Reo<  I  )  , St  1  oi  I  ) 

2330 

3  t  r  a  v  (  I  )  =  S 1 1  <  I  ;  /  S 1 1  o  (  I  ) 

2330 

NEXT  I 

2100 

CREATE  SCAT  “  STAVFCV'1  , 1  S  ,3 

2110 

ASSIGN  3Fath_5  TC  '•STAVFCV" 

2123 

1  PRINT  AVERAGE  STANTON  NUM5ERS 

c.  1  3  V 

! 

r~»  •  * 

_  , 

PRINT  "ROWS  X ' M  }  REYNOLDS  Nq  S  t / 3 1  o 

Ave 

175 


PRINT  "  Tr  LATE  -  T  FREESTREAM  CORRECTED  I  DEB  C;  =  "iDeUc 
PRINT 

PRINT  "AVERAGE  PLATE  TEMP.  CGRRECTEC  ( CEG  C  j»*" ;  Tava-Ct 
PRINT 

PRINT  “NONQIMESIONAL  TEMP.  THETA  m" t Theta 
PRINT 

PRINT  'No  X;m?  Z.;CM;  H<y/M*2  C)  St/Sto 

PRINT  "**-**»**+4**«'****'****»****'*-***-****<l***«*-*4»****4*#***< 

PRINT 

ASSIGN  SPathB  TO  “H0AT2  * 

ASSIGN  OPathT  TO  “SWC3'1 


3 1  /  3 1  f  ‘ 

f  *  *  *  4  4  ♦  « 


1  Z  1  3 

ENTER  OPathSi Ho(  I  ) 

i  :zj 

ENTER  OF  a  t  n7  5  3 1  f 1 1 

’  2  j3 

NtX  r  I 

1  240 

;CALCLLATE  THE  HEA 

•  252 

: 

;:so 

FOR  I = !  TO 

1  270 

T  ■■  I  )  =  T  ( 

4  n  ^ 

1  uZO 

H>  x  .'=0  -' 

!  290 

St-:  I  )=H 

!  202 

S  t  r (  I  )- 

i  j  » z 

3  W  r\  i.  ! 

\  3Z3 

NEXT  * 

•  JJ0 

'GENERATE  X  .AND  Z 

t  Z*'Z 

Z(  !  7 

1333 

is 

1  353 

FOR  1  =  2  TO  Z 1 

«  7"  -t 

!J-U 

X  Z  <  I  '<  =  1  .  ;  5 

i  3S3 

Z'. ; ;=zc I-i ;+! .27 

13S3 

NEXT  I 

«  *  f* 

f  *•+•*;*/ 

Zv  22  '  =  Z(  ! 

i  4  i  0 

xzcz:  :-  =  i  .25 

1423 

FOR  1-23  TO  42 

i  430 

X2C  I  >1  .25 

1  440 

Z-  I  >«Z<  1-1  >4-1.27 

1  490 

NEXT  I 

1450 

Z C  43  >Z'  1  ) 

I  473 

XZU3  >*  1  .  40 

i  430 

FOR  1-44  TO  S3 

1  450 

i  )-i  .43 

1 530 

Z '  I  )  -  Z  {  I  -  ’  1 H  1  . 2  7 

1510 

NEXT  I 

i  a  2  3 

Z  ■!  S  4  i  =  Z  (  1  1 

;53o 

X  2 ( S  4 <  * !  .5 

i  543 

FOR  1=65  TO  S4 

1553 

X 2 \  I  /  =  i  .3 

1 553 

i  >«:■:  i-i  >+i  .27 

1573 

NEXT  I 

1  5S3 

V  T  —  *  ^ 

Aj-I  .0 

*  530 

FOR  J=5  TO  5 

!  500 

n=J»2 ! -23 

iVljlj'lj’ti'tjlj 


> .» .  > .  ■ 


•jj- 


r  h  I N  •  "*■■**■*♦  +  **«*****■»■**♦*#•#*»*##►•♦♦*##***  + 

FOR  r=i  TO  5 

r".  I  N'T  US  I  NO  2  83 ;  I  ,X  i  C  I  )  ,R-y  ■;  I  )  ,  3  tra I  > 
;jiJ  ,  4/  , uu . Du  ,4X  ,u,4Cfc#7X  t  u  ,  4CE 
OUTPUT  OF  a  t  h__5 ;  X  i  (  I  5  ,  Ray ;  I  5  ,Str<av?  I  ) 

NEXT  I 

ASSIGN  3 Path  5  TC  * 


♦  ♦  ♦  *  •»  *  *  ■*  ♦  *  i 


GOTO  2560 

!  SUER OUT I ME  TO  CALCULATE  DISCHARGE  COEFFICIENTS 

! 

ASSIGN  OFatnj  TC  “FILDT'0" 

ENTER  !?Fa  i n3 ;  F i w  ,Deip  i 
1 GE i  AVERAGE  PLENUM  TEMPERATURE 
Tp  I  =3 


ipl’ i Pi/ j 


IGNVERT  R uENUM  TEMP.  TC  INJECTION  TEMP. 


O:  .50!lS. 


.  L  j  d  C  5  -  :  i  n  J  J  .  15 

U  ’  =  F  I  u/  .  30332  i  453  ‘FDR  13  INJ.  HOLES 

1  Lii  n.i=F  Iw/  .  303436  IFOR  7  IN  J.  HOLES 

1  or = ! i j ab s - i U i n j  ‘ 0 / \ 2 *  1 3GS  » ) 

Roc=PaMb/  ■'  Tcr*ZS'’ ) 

Imr  i.N=Roc*Uxn  j  !  INJECTION  HASS  FLUX 

Rop i  =  F  amb/ ! Tijabs*DS?  )  'TOTAL  DENSITY  BASED  CN  To 

L'3—  <  Da  Ip  ;  *2/Rop  I .5 

P  i*i  f  I  .v  =  R  o  g  1  *  U  3 

C  d  ■*  I  ">  f  I  /  F  m  f  i  .• 

Reb=U i n j  * , G3So2S/ . 3G03  i  EG 
M= In  f i x / (  Rc  *U inf)  !  SLOWING  RATIO 

Vrat-Uinj/Uinf  'VELOCITY  RATIO 

Mfr«Roc*Uinj “I/< Ro*Uinf “2  >  i MOMENTUM  FLUX  RATIO 

i 

!  PRINT  RESULTS 

PRINT  "INJECTION  TEMPERATURE  '  0  >«“  hip  j  PLENUM  TEMP  <  C  >-“  i  Tp  I 
PRINT 


2643 
2SS3 
2563 
2673 
25SG 
26  33 


PRINT  •’PLENUM  PRESSURE  DIFF,  •:  N/M'2  >="  sOalci 
PRINT 

PRINT  ■’ COOLANT  OENSITY  KG/M'3  ;  =  “;Roc,“  DENSITY  RATIO  =  "  ;Roc/Rc 
F  R I N  T 

PRINT  "INJECTION  VELOCITY  (MY  3  )  =  “  iUin  i  ."MASS  FLUX  <KG/M"2  S  I  m  f  1  > 
PRINT 

Pr.  i  NT  Rc.:iJuLuo  No  '  Cin  ia  jRab.^CiaCHARGE  CuEFF.  Ccj  = "  ;  C  d 
PRINT 

PRINT  "SLOWING  RATIO  -"jM,-  VELOCITY  RATIO  =":Vrat 
FRINT 

DRINT  " MOMENTUM  FLUX  RATIO  =  "SM?r 

PRINT 

GOTO  533 

MASS  STORAGE  IS  “  :  INTERNAL  ,4. " 

END 


% 

■y^yyy'^-' 


REM  PROGRAM  nCQTrnu 

REM  THIS  PROGRAM  ACQUIRES  TEMPERATURES  FROM  A  THERMOCCOUPLE 
REM  M.ONTED  ON  A  AUTOMATIC  TRAVERSE  MECHANISIM  AND  FROM  A  FREE 
'  i  STREAM  VELOCITY  THERMCGPLE  TO  PLOT  THE  TEMPERATURES  PROFILE': 
!  USING  PROGRAM  FLQTDE 

t 

i  9Y  ALFREDO  ORTIZ  .  AUGUST ,1987 
REM 

REM  VARIABLE  NAMES 
REM 

REM  E\ I  ■  IS  THE  VOLTAGE  REAC  FROM  THE  DATA  ACQUI3TIQN  SYSTEM 
REM  Till  IS  THE  CONVERSION  FROM  VOLTAGE  TO  TEMPERATURES 
REM 

-EM  =AMc=AMSIENT  PRESSURE 

REM  TFS-rREE  STREAM  TEMP: KELVIN < 

REM  RQ*CENSI7Y  vKG/M''3) 

REM  U INF  =  VELCC I T'i  Or  THE  rREESTREAN 


1  :2  c- 

DIM  v: 

,  ij  is  .  , 

Z { ’ 222 ) ,T( 1 222  )  . T 1 i ! 323 

’<  ,D<  ! 

x.  o  c 

REM 

7  «  ^ 

REM 

4.  —  C 

T  ^p;  ■  J 

’  E  n  7  E  ~ 

DATE  ,  (  MMOuYY  } "  ,fcS 

AMA 
c  -  V 

INPUT 

" ENTER 

TIME  ,(KHMM  V  ,N3 

Z40 

REM 

zso 

IMFLT 

”  PNTEP 

DISTANCE  FROM  THE  3.L. 

TR  I  r 

ZS3 

INPUT 

"  ENTER 

POINTS  5PANWISE ”  ,M3 

272 

INPUT 

“ENTER 

POINTS  VERTICAL  (MUST 

3E  AN 

zsc 

REM  N3 

MUST 

SE  AN  EVEN  INTEGER 

ZSG 

INFUT 

“ E  NTS R 

SFANWISE  RESOLUTIONS  IN 

)  ,  *.  4- 

322 

I  Nr  LT 

“ ENTER 

VERTICAL  RESOLUTIONS  IN 

/"  ,Y4 

Z !  3 

INPUT 

•'  I N I T I 

AL  I ( IN  T  ,Z2 

A 

*-  u 

INPUT 

" INITI 

AL  Y\ IN  ;  “  ,  Y3 

REM 

Z( !  )  =  Z3 
y>  i  )«y: 
N7=N3/  Z 


J  I  =Z*IS’*-M3-3 
J  Z  =  Z  *  I S  *M3rH3 
FOR  K - J i  ID  J Z 
Z(K >-Z< « - >  )+Z4 

Y  (  K  i  *  Y  (  K  -  i  ) 

NEXT  K 

JD-Z* IS  *MI  +  M3+ 1 
Z !  J  3  ;  =  Z J  3-  i  5 
Y '  J  3  /  •  Y  S  j  3  -  ’  'r  Y  - 
J4=J3t! 

J5=2  » IS  *-M3rZ  *M3 
FOR  K  =  J  4  "0  J  5 
Z  \  K  )  =  Z(K-I  --Z4 

Y  K  5  =  Y K  - !  ) 

NEXT  K 


jwVjv&.V.v'Ov'CvVV _%*_««  •.•-jv' v. 


540 

550 

560 

570 

530 

530 

600 

6  i  0 
GZ0 
630 
640 
650 
650 
S70 
630 
5  33 
73C 

7  i  0 
720 
73S 
743 
750 
762 
773 

•->  r\  r* 

i  ov 

730 

600 

813 

323 

SZ2 

640 

353 

363 

3r- i  r\ 

:  w 
330 
390 
900 
910 
923 
930 
940 
950 
950 
970 
930 
990 
1000 
1010 
1020 
1030 
1043 
1050 
1 060 
1070 


IF  I7*N7  THEN  590 

J  6  =*  J  5  + 1 

Z<  JS  >  =  Z(  J6-1  ) 

Y  (  J  6  )  =  Y  (  J  6  - 1  / + Y  4 

NEXT  17 

REM 

i  ENTER  FREE  STREAM  TEMPERATURE  PROM  DAS. 

OUTPUT  709 ; " A I " ; 151 i’UTI' 

ENTER  709; X 
E-X* 1000000 

7f  5=-.  0733  171  •<■.3261 0S4S4«»E-.  000000678  *E*E 
T  f 5  =  1 f s+273. 1 5 

REM  ENTER  AMBIENT  CONDITIONS 

INFUT  " ENTER  FAM3  (IN. OF  KG;"  .Fane 

INPUT  “ENTER  PP.E33URE  OIFFcFENCE  (IN  CF  H20  .Da:,  tap 

REM 

REM  CONVERSION  TO  31  UNITS 
Pamb=F5!yib*333S .  32 
Ro**Pai*ib / ( 237*T  f  3 / 

Deitap=0eltap*243.7 
REM  FPEE3TREAM  VELOCITY 
Ui  n  f  =  (  2  *DeI  t  ap/Ro  •'  ‘ .  5 


REM 

REM  ENTER  THE  LOOP  FGR  ACQUIRING  EACH  TEMFEF.A  ,•  URE  COMFU 
REM  AND  CONFUTING  TWO  TEMPERATURE  COEFFICIENTS  AND  TUG 
PRINTER  13  701 

PRINT  CHR3CZ7  )jCHR3< 36  >;CHR$< 106  )?CHR$( 53  );CHR3( 52  )  ;CHR 
PRINT  CHR3 ( 27  > ;CHR3( 33  )  ;CHRS( 108 ) ;CHRS( 49  )  jCHRSC 7c ) 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


“TEMPERATURE  PROFILE  COMPUTATION  " 

"DATE  OF  RUN  IS“,N8,”  TIME  =";N9 

“DISTANCE  FROM  THE  BOUNDARY  LAYER  TRIP.  X  (  M  )  =  " ; X 

“DENSITY!  KG/M-'p  5“  ,fic 

" FREE3TREAM  VcLOCIiY  (M/S/“  .Uinf 

“  PAMS  I  ENT  t  N/M  “  2  / "  ,  Paf'io 

"FREE  STREAM  TEMPERATURE  (DEG  Ci  =“ ; Tf s-273 . 1 5 


PRINT 

PRINT  “POINTS  3FANWISE" ,M3 
PRINT 

PRINT  “POINTS  VERTICAL" ,N3 
PRINT 

PRINT  “SPANWrSE  RESOLUTION  (IN)“,Z4 
PRINT 

PRINT  “VERTICAL  RESOLUTION  IN)“,Y4 
PRINT 

PRINT  “INITIAL  Z( IN  ,Z3  ; 

PRINT 
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'.AT.’J'.  VA".VA^JVW'iVkV.'-'  ."V  V>'  •.'*•** , 


■jror=jr»jr'y=>7jr">  >  r  s  ■  v  .  v.  v.v>-. 


!O30 

PRINT  '“INITIAL  Y<  INK  ,Y3 

1030 

FRINT 

1  100 

PRINT  ”*  +  '»**-♦**«**•*■**♦  +  '♦**''*****♦****♦•****•♦  +  **'* *  *  * " 

I  110 

PRINT 

1  120 

PRINT  Y  Z  T  TF3 

DT" 

1  !  30 

PRINT  “=====»=*=3— ==»»«*==— =“=-*=**=—=»====•« 

sa aeaasas  “ 

1  U0 

PRINT 

1  150 

DI5F  “'HIT  <( CONTINUE  FOR  OAS>>" 

1  150 

PAUSE 

1  !  70 

K3=0 

1130 

K2=M3*N3 

1  130 

FOR  K  = 1  TO  K2 

1 200 

K3*K 

1210 

REN 

1223 

WAIT  13 

•  —  ~J  w 

REV1  ACQUIRE  THE  TEMPERATURE  FOR  EACH  POSIT 

I  0»v 

t  2.  •+*/ 

61=3 

1.253 

REM 

1  2  S3 

FOR  J  =  1  TO  25  !  ENTER  THE  GAS  AND  SAMPLE 

OF  TEMPERATURE 

1273 

OUTPUT  703  ;  “ft  I” 152  i  '"v'TI  " 

1230 

ENTER  703  iX 

1  233 

G I  *5  !  +  .< 

«  “!  -S  f\ 

1  J  1/  L 

NEXT  J 

1313 

REM 

1  32C  X=G1/2S  ! AVERAGE  THt  VALUES 

f  3ZC  E  =  X  *  1 333230 

t  343  T( K  )=.  1  S3S  733+.  325557!  3  *E-.  080302488 -*fc*E  ! CONVcRS ION  FROM  VOLTAGc  TC 
PERA7URE 


1 353 

REM 

!  352 

1  Sz.~~  FREE  STR 

1  370 

G4=G 

1380 

FOR  I = 1  TO  5 

1330 

OUTPUT  73S;"AI 

U00 

ENTER  739 jX 

1410 

i34=G4rX 

1420 

NEXT  I 

1  430 

X=G4/5 

1443 

E=X * ' 333000 

I  453 

T 1 ( K 1=- , 375 j 1 7 

1450 

D  '1  K  /  *  T  \  K  1  -  T 1  (  K 

1  470 

PRINT  USING  14 

1430 

IMAGE  MOD . DO  ,2 

1  430 

REM 

1530 

NEXT  X 

1510 

REM 

1  520 

MASS  STORAGE  I 

1  530 

CREATE  BOAT  “ T 

1  540 

ASSIGN  ?Path2 

1553 

FOR  1*1  TO  K9 

*  550 

OUTPUT  SFathZ; 

1570 

NEXT  I 

1530 

ASSIGN  OPatnZ 

1530 

MASS  STORAGE  j 

1 530 

END 

JLtSSSJSJ. 


I  PROGRAM  FLOTPRUN 

!  THIS  PROGRAM  READS  THE  FILE  TPRO  AND  PLOT  TEMPERATURES 
t  FRCFILES 

DIM  /I  C  300  )  , 23(  500  ;  ,XS<  800  >C  II  , LM  (  800  )  ,ZI  {  40  ,*  ,Z2<  40  ;  ,C3(  13  i[  1  3 
!  U(I>=D(I>  :  TEilFERATURE  DIFFERENCE. 


INPUT  "ENTER  RUN  NUMBER  ''  ,Rufino 
FS  = " DATE= " 

MASS  STORAGE  IS  ” : INTERNAL  ,4 , 1“ 
ASSIGN  @Pathl  TO  “ TFR07 " 

FOR  I-i  TO  300 


ENTER 

@?at  h 1  : V 1 ( I  1 

,201 I  1  ,U1( 

I  ) 

V  }{  l 

=7!  (  I  )•*"•.  S4 

IG‘  I  * 

=zo! :  >*:.54 

MEX  T 

i 

eQR  I 

-1  TO  300 

if  u; 

s  I  THEN 

X  3  •;  I  '  =  ■'  3 " 

IF  \J 1 

1  •:  I  ).;  .2  TH 

EN  X3 

:  I  >  = 

"  1  •• 

;  F  V  - 

il  ) <  =  1  .3  AND 

u  i  ■:  i  »  .  5  : 

HEN  X3< I  > 

a  '*  2  " 

IF  U! 

i  ,5  AND 

v  1  v  I  1 /  i  .2 

THEN 

X  3  ( I 

I  ~  V  • 

;lK=:.0  AND 

V  ’  •  1 1 ;  {  .  5 

THEN 

XS4  I 

if  u; 

( I v  \  =  2 . 5  nNC 

V1 !  '  I  ;  ✓  Z  .  0 

THEN 

X  5  {  I 

)  = "  5 " 

if  u: 

/  -r  •  T  '  r*> 

•.  i  ;  .-  j  .  'a  r\*u 

V 1 V I  >  :•  2 . 5 

then 

x  5  •:  I 

)  -  ”  c  N 

I F  v  i 

i  I  ) <  = 3 . 5  AND 

'Jill )  /  j  .  0 

them 

X3-:  i 

i  ~  i 

IF  U i 

(IK-4.0  AND 

U1  (  I  >33.5 

then 

X3<:  I 

;*"3“ 

IF  U! 

i  I  .'<  =  4.3  AND 

'j  i  ( I ;  >  4 . 3 

THEN 

XS(  I 

)  =  “  3 M 

IF  v! 

;i><=5.0  AND 

') }  {  I  1/4.5 

THEN 

A  S  ;  I 

! —  a  ’ 

IF  U1 

f  I  ) '.  =  5 . 5  AND 

OKI  )  ■  5 . 3 

THEN 

x  s  ■;  i 

)  *  "  b  “ 

IF  V! 

( I  1/5.5  THEN 

X  3  (  I  )  = "  c  “ 

NEXT  I 
3 1  NIT 

PLOTTER  IS  70S  ,  "HFGL " 
i PLOTTER  IS  CRT  ,  "INTERNAL" 

GRAFHICS  ON 
C3IZE  2.5, .35 
MOLE  35  ,  i A 

LABEL  "EMBEDDED  LOR TEX  10  M/S  WITH  FILM  COOLING,  M=3.98" 
MOVE  35 , 1 1 

L.-.EEl  '  a  —  1  .  43G  M  ,  UNHEA  i  EC  PLATE  ,UGRi EX  GEN.  AT  4.73  c.^> 
MOLE  75,33 

LABEL  USING  "SA  , 4  ,60 . 40” ;F3 ,Runno 
031 ZE  4 . 5  , . BE 
MOUE  40  ,  So 

LABEL  "TEMPERATURE  PROFILES" 

CSIZE  3 .5  .  .65 
MOUE  57,17 
LABEL  Z  <  CM  >  * 

DEG 

LOIR  SO 


!  PROGRAM  PTSTAV 

! THIS  PROGRAM  READS  A  DATA  FILE  AND  PLOTS  AVERAGE  STANTON  NUMBERS 
l VERSUS  REYNOLDS  NUMBERS. 

i 

DIM  Rey (  5  ;  ,X1<S  )  .St:  (S)  .CSC SI  ,St<S  ; 

MASS  STORAGE  IS  *:  INTERNAL  ,4. .  1‘ 

ASSIGN  §Path2  TQ  “STAVI " 

FOR  I = !  TO  S 

ENTER  SPathE ;  X 1  \  I  :  ,Rey<  D.StIC) 

NEXT  I 

PRINTER  IS  I 

PRINT  "ENTER  RUN  NUMBER  TO  WHICH  'HIS  DATA  CORRESPOND  >  MMODV.'  . HAMM  > 
INPUT  Runno 
F3=  ’DATE  «" 

GINIT 

PLOTTER  IS  705  , “HPGL" 

! PLOTTER  IS  CRT  ,  “  INTERNAL" 

GRAPHICS  ON 
CSI2E  2 .5  ,  .65 
MOVE  40  ,  U 

LABEL  “ FREESTREAM  IS  M/S,  NC  FILM  COOLING" 

CS 1 2E  4 . 5  .  . 55 
MOVE  43,35 

LABEL  AVERAGED  STANTON  NUMBERS" 

CSIZE  0.5,  .55 
MOVE  75,53 

LABEL  USING  "5A  ,SD . 40" ;FS  .Rubric 

CSIZE  3.1..S5 

MOVE  6  2  , !  7 

LABEL  "RE*  *  ! 3 " 5 " 

DEG 

LOIR  30 
MOVE  20,33 

LABEi-  "STANTON  No  *!0"-3" 

LOIR  3 

VIEWPORT  30,  MS,  OS, 33 
WINDOW  0  ,5  ,3  ,S 
FRAME 

AXES  .5,.  5, 0,0  ,2,1 
AXES  .5  ,  .5  ,5  .5  ,2  ,1 
CLIP  OFF 
CSIZE  0 .3  ,  .65 
LQRG  2 

FOR  1=0  TO  S  STEF  .5 
MOVE  -.22,1 

LABEL  USING  "#  ,0.00" ;  I 
NEXT  I 
L  0  F.  3  4 

FOR  .J  =  0 .  TO  5.0  STEP  .5 
MOVE  J  ,-.22 

LABEL  USING  "#  ,00.0"  ;  J 
NEXT  J 


v"<’vvv',,>r.Tv".\,.’V’^;<".V-V'.'^ -\x\»  wu«  r*r>  =j.-iar»ir><  •.-*  -■>  .vv 


530 

'ENTER  VALUES 

FOR  ANALVTI 

CAL 

FORMULA 

54.0 

Zi-I . 1 

'METERS 

550 

Pr=. 71 

1 PRANOTL 

No 

FOR  AIR 

550  ‘ 

FOR  1=1  TO  6 

570 

St( I )=( ( .030* 

(  Rey  ( I  . 

2  )  ) 

*<  1  Zi/Xl < I  )  )A 

5BQ 

S.t(  I  )=5t (  I')*t 

.  0E+3 

5-30 

Rey< I  )°Ray\ I ) 

♦t . 0E-S 

500 

PLOT  Rev < I  )  ,S 

t  <  I  -* 

S.10 

NEXT  I 

520 

CSIZE  1 .3, .72 

530 

FOR  IH  TO  5 

54.0 

C$CI]="x" 

550 

st  m  i  ,=st i : i ) 

*  1  .  0E*Z 

553 

MOVE  Rey  •;  I  ;  ,S 

t  !  (  I  ) 

570 

LABEL  USING  " 

A "  i  CSC  I  ] 

530 

PENUP 

533 

NEXT  I 

700 

MASS  STORAGE 

IS  INTERNAL  .4 

« 

713 

END 

184 


'luiiiii.iliimi. 


.V^^./ClvVa.V. 


10  ! PROGRAM  PTSTLC 

20  'THIS  PROGRAM  READS  A  DATA  FILS  AND  PLOTS  SPANWISE  VARIATIONS 
30  'OF  HEAT  TRANS  PER  COEFFICIENTS 

40  i 

50  DIM  H<  200  )  ,St<  200  )  ,Rey(S  )  ,21 230  )  .XI  <6  >  ,X2(200  >  ,B(  1  25  )  ,C$<  c  it  3 

60  MASS  STCRAGE  IS  *: INTERNAL  ,4 , 1 " 

70  ASSIGN  3Path!  TO  “H0AT2" 

80  FOR  I-l  TO  1 2S 

90  ENTER  SPathi :H( I >  ,St( I  )  ,X2< I ) ,Z( I ) 

100  NEXT  I 

110  F$="GATE  ='• 

123  INPUT  "ENTER  RUN  NUMBER  ( MMCCYV . HHMM ) " .fiunno 
’23  3INIT 

140  PLOTTER  IS  735  /HPSL" 

150  'PLOTTER  IS  CRT  , " INTERNAL ” 

ISO  GRAPHICS  ON 

170  CSIZE  2.5  ,  .65 

130  MODE  42,14 

133  LABEL  "FREE  STREAM  10  M/S ,  NO  FILM  COOLING' 

203  CSIZE  4.5,  .65 

210  MOVE  22,65 

223  LABEL  "  3PANUI I SE  HEAT  TRANSFER  COEFFICIENT- 
230  CSIZE  2 . 5  . .  S  S 

240  MOVE  75,32 

252  LABEL  USING  "6A  ,4  ,60.40"  ;F3  , RLir.ro 

2S2  CSIZE  2. 5, .65 

270  MOVE  55 , 17 

230  LABE'-  "Z  (CM)  * 

230  MOVE  22,25 

300  OES 

310  LOIR  93 

320  LABEL  “H  U/M'2  DEG  C* 

330  LOIR  0 

340  VIEWPORT  30 , 1  IS  ,25  ,83 

350  WINQCU  -15.24,15.24,0,100 

3B0  FRAME 

373  AXES  1 .27 , 10  ,-15.24  ,0  ,2  ,2 

380  AXES  1.27,13,15.24,100,2,2 

390  CSIZE  2.0,  .65 

400  MOVE  -13,23 

410  LABEL  "X  ( M  > " 

420  A=22 

430  FOR  K= 1  TG  6 

440  D  =  A-r;*3 

450  CSC  1  i*"-X=1 . 15" 

460  C$( 2 >« "=X»I .25" 

470  CSC  3  )="=:<- 1  .40" 

430  C$C4)="=X*!.S0" 

490  CS(  5  )=•  "=:<=  1  .  93 " 

500  C3cS  )-■'-/=>  2 . 00  ” 

510  MOVE  -14  ,0 

520  LABEL  USING  "0  ,8A" ;K  ,CS\  K  it  1  ;3  ] 

532  PENUP 

540  NEXT  K 
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CSIZE  1 .3  ,  .72 
FOR  J » !  TO  S 
M-J*21-20 
FOR  I=M  TO  J*2 1 
B< I  )  =  J 

MOUE  Z< I)  ,H( I  > 

LABEL  USING-  "O'  ;B<  I  ) 

PENUP 
NEXT  I 
NEXT  r 
CLIP  OFF 
CSIZE  2.0..B5 
LQRG  2 

FOR  1=0  TO  103  STEP  10 
.  MOUE  -17.75,1 
LABEL  USING  "t  ,DDQ . Q" ;  I 
NEXT  I 
LORG  S 

FOR  J=- 1 S . 24  TO  >5. "4  STEP  5.08 
MOUE  J  ,-2.3 

LABEL  USING  "4  , ODD. DO" ;  J 
NEXT  J 

MASS  STORAGE  IS  " : INTERNAL  ,4" 


r RuaRnrl  rlo  iU 

THIS  PROGRAM  READS  A  QATn  FILE  AND  PLOTS  THE  RATIO  OF  STANTON  NUMBER: 
OV/ER  BASE  LINE  STANTON  NUMBERS  UERSuS  REYNOLDS  NUMBERS . 

IK  Rev i  =  /  ,X  i  ( S  )  ,S 1 1  ( 3  ) ,  Si Z 1 5  >  ,Re  1  (  S  >  ,C${  o  )£  1  I  ,?S‘  E  ) I  !  j  fGS<  a  ; l  1  J 
I M  S 1 3  C S  )  ,  ReZ ( S  /  ,  S 1 4 \ S / , R e  4 C S / , 3 1  r I <5 ) , S  t  rZ < B I , S  t  r 3 '  G ) 

SSISN  ?Fathl  TO  "STAv' 

GR  I a I  TO  S 

N7ER  3Pa  t  h I  ;XI i I  >  ,Fey  <  I  /  ,S  t i ( I  / 

EXT  I 

33I3N  QFatrrZ  TO  "STnUFC" 

OR  I = ’  TO  S 

NTER  3Path2?Rei<  I  )  ,3i2-:  I  ) 

tr  !  ■!  I  / =  3 1 Z 1  i  //’Sti'  i  ) 

S S 13 N  ?Fatrij  TO  "bi'iUFl 
OR  I*!  TO  5 

NTER  3Fa  t  n3 :  Re3C  I  '  ,  3 1 3 ■'  I  ) 
t rZ I  :  =5i3  ;  I  )  /St*  i  I  ) 

:e:xt  I 

33ISN  ;0Patn4  TO  ”STAUr2" 

OR  !=■'  TO  S 

NTER  5»ath4;Rs4t  I  >,SU(  I  > 
t  r  3 ;  I  =  3 1 4 I  >  >  S  t-  i  •  '  I  ! 

;=;<T  I 

■R  INTER  13  1 

PINT  "ENTER  TODAY’S  CATE  ( HMOOVY .HHMM  ;  ' 

INPUT  Run  no 
5-  "DA  •.  E  = 

HINT  'ENTER  IN  ORDER  CORRESPOND  IMG  BLOWING  RATIOS,  Ml  ,KZ  ,  AND  M3” 

OF.  I  —  i  >C  j 
NPUT  ri  <  I  ) 

ECHO  PRINT  M 
PINT  ”M»" ;M< I > 

UNIT 

PLOTTER  IS  7eS.-HP0L* 

'LOTTER  13  CRT, -INTERNAL  * 
iPAFHICS  ON' 

‘_30u_i'ias  *  t  GO *MAX,C  !  ,HATIO 
r_gdu_na.-: “  I  00*MAY  11  ,  1  /  RATIO  i 
IS  I ZE  4  .  S  , . S5 
.OPS  S 

IDLE  X_gdui_r,ia.'. /Z  . .  9’ 

.ABEL  “FILM  COOLING"  - 
5E3 

.DIP  30 
:si:e  3.b 

10UE  .  !  *Y_gdu_na.'.  ,Y_gdu_m<4*/2 
.ABEL  "St/Sto" 

.DIR  3 
.ORG  4 

1GUE  .<_gdw.  _roa.-. /2  ,  .  Z0-* V_gdu_Ma.< 

.ABEL  ’RE /MOLDS  No*!3’’S" 


V  "/  /.V.  «' 


•  ‘.'VvVsT 


s’.v'.Y.s'.sV 


'PROGRAM  PL3TRV0R 

(THIS  PROGRAM  READS  A  OATA  FILE  AND  PLOTS  SPANWISE  VARIATIONS 
'OF  STANTON  NUMBERS  RATIOS  BY  ROUS 


DIM  Str( 233  )  .Rey! 6 ) ,Z( Z0O ) ,X1 <S  )  ,X2< 200  > ,CS( S  >C8 1  ,8< 1 25  > 

MASS  STORAGE  IS  * :  INTERNAL  ,4,1" 

ASSIGN  ©Path !  TO  ‘ STRDT2 “ 

FOR  1=1  TO  1 26 

ENTER  ©Path ! «5tr( I ) ,X2(  I > . Z < I > 

NEXT  I 
F$*"OATE  ='• 

PRINTER  13  1 

PRINT  "ENTER  RUN  NUMBER  TO  WHICH  THIS  PLOT  CORRESPOND  < MMDDYY 
INPUT  Runno 
6 1 N I T 

PLOTTER  IS  795  ,  "HFGL “ 

(PLOTTER  IS  CRT  ,  ” INTERNAL" 

GRAPHICS  ON 

.’C’-c  i  r  e  — 

MOVE  42.14 

LABEL  "FREE  STREAM  10  M/S,  NO  FILM  COOLING  " 

CSIZE  4 . 5  , . 55 
MOVE  23,35 

LABEL  "STANTON  NUMBER  RATIOS,  VORTEX  GEN  52" 

C3I2E  2. 5, .55 
MOVE  70,33 

LABEL  USING  h SA  ,5  ,50 . 4D"  <Fi  .P.unr.o 
CSIZE  3.5  ,  .55 
MOVE  55,1? 

LABEL  "Z  <  CM  >  " 

MGVE  22  ,45 
DEG 

LDIR  90 


LABEL  "St/Sto" 

LDIR  0 

VIEWPORT  30,115,25,83 
WINDOW  23  ,-20  ,  .4  ,! .4 


FRAME 

AXES  5 , . 1 , 

23,. 

AXES  5  , .  1  , 

-20  , 

CSIZE  2.3, 

.55 

MOVE  13.5, 

.85 

LAeEL  "X  i. 

Ml" 

A= .  83 

FOR  K= 1  TO 

6 

D=A-K ♦ .04 

CS(  1  )= " =X  = 

1.15 

C$i 2  )  =  " =X  = 

1  -25 

C  3  (  3  )  = "  =  X  = 

1  . 43 

CS( 4  )  =  " =X  = 

1  .53 

C3' 5  /  = " =  X  = 

1.80 

C  £  (  b  '  =  "  *  X  = 

2.00 

MOVE  14,0 

LABcL  USING  *0  ,3 A" ;K  ,C3(K ) 

PENJP 

NEXT  K 

CSIZE  1.3, .72 
FOR  J=l  TO  S 
M=J*21 -20 
FOR  I»M  TO  J *2 1 
B< I  )«J 

MOUE  2( I ) ,Str< I  ) 

LABEL  USIN6  “0“ ;B( I  ) 

PENUP 
NEXT  I 
NEXT  j 
CLIP  OFF 
CSIZE  2.0, .55 
LORG  2 

FOR  I=. 4  TO  ;  .4  STEP  .  I 
MOUE  21.7,1 
LABEL  USING  '4,D.C‘:I 
NEXT  I 
LGRG  5 

FOR  J=-20  TO  20  STEP  5 
MOVE  J  , . 37 

LABEL  USING  " z  ,000 .00"; J 
NEXT  J 

MASS  STORAGE  IS  INTERNAL  ,4 


1  PROGRAM  PlSTRWS 

'THIS  PROGRAM  REAOS  A  DATA  FILE  AND  PLOTS  SPANLISE  VARIATIONS 
! OF  STANTON  NUMEERS  RATIOS  BY  ROWS  FOR  FILM  COOLING  DATA  AND 
'FOR  VORTEX  OATA 


DIM  Strl <  )2G)  ,St r 

2(125)  ,Z(  IZS  )  ,y.2 

D I M  S  t  r  t  i  2  5  '  ,  2  i  ;  1 
| 

2S  ) 

MASS  STORAGE  IS  ” 

: INTERNAL  ,4,1” 

ASSIGN  2F3 t h 1  TO 

" STRFC3" 

ASSIGN  ;?patn2  TO 

” STRFOV 1  :  ” 

FCF  !«;  TO  ’05 
EN'ER  0“3 : A ' ; G t  —  ’ 

■  :I  i  ,X2(  I'i.Zil  > 

ENTER  Tv  a  t  r\Z  :  5  trZ 

;  i )  ,«:?•  i :  ,xz; :  > 

ri=  •DATE  =” 

PRINTER  IS  i 

FR  IN"  ” ENTER  FUN 

NUMBER  TO  WHICH 

a  I':  'J  ,  r  , ;  n  D  C 


i  PLOTTER  IS  CRT  , ” INTESNAu " 

GRAPHICS  ON 
C3IZE  2.5  ,  .55 
MOVE  :o,u 

LAEEL  ” -REE  STREAM  15  M/S  ,  WITH  FILM  COOLING,  VORTEX  GEN.  AT  4.75  cn 

Q  C.  I  ?  £  s  C  _  C  £ 

MOVE  40,05 

i_AEE_  -  STANTON  NUMBER  RATIOS” 

CSIZE  :.5,.SS 
MOVE  75,33 

lABEl  USING  "3 A  ,Z  .SC.AC  ‘  :FS  .Eunr.c 
CSIZE  Z.5,.55 
MOVE  S 5,1  7 
LABEL  Z  ; CM ;  ' 

MOVE  05,45 
DEG 

LOIR  90 

LAEEu  'St /Etc" 
lDiR  j: 

VIEWFCR  :  j  0  , :  *  E  .  2  s  ,  S  _*■ 

WINDOW  -00  ,20  ,  .3  , : .2 

::  a-^c 

,->ES  5  ,  .  !  ,-20  ,.30,!  ,! 

AXES  5  . .  i  .20 .1.2.1  . ! 


A0-A193  1SS 


THE  THERMAL  IEHMVIM  OF  FILM  COOLED  TURRULENT  ROUNDARV 
LAVERS  AS  AFFECTED  RV  LONGITUDINAL  VORTICEStU)  NAVAL 
POSTGRADUATE  SCHOOL  HONTEREV  CA  A  ORTIZ  SEP  A? 

F/G  21/13 


UNCLASSIFIED 


! PROGRAM  SURFCONT 1 

! THIS  PROGRAM  READS  A  DATA  FILE  AND  PLOTS  SURFACE  CONTOURS  OF 
!CF  STANTON  NUMBERS  RATIOS  ALONG  THE  PLATE  SPAN. 


i ORTIZ  VERSION  JULY/ST 

DIM  St-(  5  25  >  ,Strf  U2S  >  sZv  126  ,X2<  I2E  :•  ,B<  iCS  >  ,Zi  t  5  '  ,I2l5) 
MASS  STORAGE  IS  *: INTERNAL  ,4  ,  I " 

ASSIGN  ©Path!  TO  “STRFCV12" 

PRINTER  IS  70! 

FOR  I-i  70  ! 25 

ENTER  ©Pa: A !  s Str(  I  •  ,3tr?C I '  ,X2<  I  ,Zi  I  ? 

NEXT  I 


lie 

SCT.VT^P  T  C  1 

;  cn 

PRINT  enter  run  Nur 

EEF  TG  WHICH  THI 

1  GO 

INPUT  Run re 

170 

» 

!  50 

5  t  r.ir*  .  5? 

1  SC 

i 1 ;r* . 55 

—  w 

FOR  !=■  re  : :s 

*“  •  .  > 

1  V.' 

IF  Str  f  <.  I  :■  -.Li  ir*.  THE! 

6v  I  =0 

I-  Sr-f  •  I  '•■>!  .32  ANZ 

Str r ! I  ;>L i i~  In 

AID 

IF  S t r r’ ■;  I  >•'=■;.* C  A?iE 

3 1  r  f  i  .  C  _  id 

4!  i 

IF  Str? <  I  *j  .20  AND 

3 1  r  f ,  I  /  ."•  i  .  !0  TH 

AID 

IP  Str?-:  I  ?>  !  .ID  THE' 

S\  I  )  =  4 

*  T  ,*i 

NEXT  I 

4  50 

SIN  IT 

453 

r  lOtT5! 

470 

i  F’l.0T7! 

A3C 

GRAPH  I  • 

450 

C  S I Z  £ 

300 

MOVE  ! i 

510 

-ABE. 

EZ0 

CGIZE 

S3  j 

MOVE  4 

=  40 

LABEL 

5=0 

CSIZE 

5b  0 

MOVE  7 

£73 

LABEL 

ESC 

CSIZE 

EGO 

MOVE  5 

SCO 

-AEEL 

£  !  3 

MOVE  i 

=  Tit, 

0E3 

£33 

uDIR  S' 

o  4C 

LABEL 

650 

LOIR  0 

Sc  3 

CSIZE 

£73 

MOVE  S' 

£83 

LABEL 

zzi'i .  A I 


’  x  •:  r  t )  n 


i  JV  A  A  -.V  -N  '  .  *»  A  '  .  »  .  »  m  «  ■  *  a  »  ■*  i  « 


«J  «• 


sso 

A=24 

730 

FOR  K =0  70 

710 

3=>A-K*2 .5 

720 

MOUE  34 ,5 

730 

Z2i  0  )*=5t.vun 

740 

Z !  <  0 .' »  3 1  n  i  n 

750 

Z2i  !  .‘=.58 

753 

ZU  !  )=!  .  32 

773 

z .3: 

750 

z:  2  •= !  .  ’ o 

7SQ 

LABEL  USING 

30G 

I0A5E  D  2'- 

i  ,-v 

bZvy 

-or-  Tu 

SJG 

3-r.-1  <~Z  ‘  *2 

r.  Jv; i  o 


LAS5L  US 1 NS  SS0:F  ,22'.  !•,  i  , 2  M  K  > 
i\£-  '  K 

u I NC-OU  ! . 1 O  ,2 .20  ,23.  ,-20. 


333 

34  2 

hSfi  V  ’-2 

a  c  3 

.35  ,5. J . 

,13  22. 

,4  .  1 

350 

SE3 

/iXEE 

nCvE 

•35.5. ,2. 
1  .  30 , 17 

.22  ,’-20. 

■  .4,1 

* 

L*6EL  -M  =0.56' 


383 

'  2  T  ^  Z  j  "7 

rn 

■  c. 

330 

FOR  I = !  70 

126 

1000 

mole  i  • , 

Zi  I  > 

1013 

LABEL  US  INS 

"0 ■ s5<  i ; 

1023 

PENL'F 

1030 

NEXT  I 

!  040 

CLIP  OFF 

1  253 

3  S I Z  -  i  3 

65 

i  3E0 

LORO  C 

1073 

FOR  I =-23  7 

0  20  STEP  3 

i  353 

MObE  !.0i 

,  r 

1030 

LABEL  US  I 

NS  “ 4 ,DCD .0"  ;  I 

I  !  03 

.  y  t 
i 

!  !  :3 

LG  PC-  3 

l  1  20 

FOR  J= i  .10 

7C  2.30  STEP  .  1 

!  1 33 

MOVE  J  ,2 

2 

i  1  43 

LAE El  US 

IK'S  "4  , 0 .  OD "  i  J 

1  150 

NEaT  j 

1  !  53 

AA S3  STORAGE  IS  * : INTERNAL 

1  172 

ES^. 

ju!d 


■■./  V.  «" 


gss«gwsg^y>y 


!  PROGRAM  FLSTRFC 

ITHIS  PROGRAM  REAQS  A  DATA  FILE  AND  PLOTS  SPANUISE  VARIATIONS 
i OF  STANTON  NUMEERS  RATIOS  3V  ROUS  FOR  FILM  COOLING  DATA 
! 

DIM  Sir <230  )  ,Rey ( 6  )  ,Z(  200  )  ,X  1  <  5  )  ,X2(  200)  ,CS<S  >£8I  ,8<  t  ZB  ) 

MASS  STORAGE  IS  ' : INTERNAL  ,4,1“ 

ASSIGN  QPathl  TO  “STRFC" 

FOP.  I- 1  TO  1 25 

ENTER  SPathl ;Str< I  )  :X2( I  )  ,Z( I > 

NEXT  I 
FS“" DATE 
PRINTER  IS  I 

PRINT  "ENTER  RON  NUMSER  TO  IJHICH  THIS  PLOT  CORRESPOND  (MKCDW.HHMM 

INFL”  Funno 

sin:: 

PLOTTER  IS  ?CE,‘HFSl‘ 

•  plotter  is  cft  .-internal- 

GRAPHICS  ON 
MOLE  -1C,  14 

LABEL  "FREE  STREAM  15  M/S,  WITH  FILM  COOLING  " 

CSIZE  4 .5  ,  .55 
MOLE  40 . 3S 

l.AEEl  "STANTON  NUMBER  RATIOS  ' 

CSIZE  :.S,.SE 
MOVE  75  ,33 


-70 

LABEL  USING  “6.A 

280 

CSIZE  3 . 5  , . SE 

220 

MOVE  55,17 

300 

LAEEi-  “  Z  (CM) 

310 

MOVE  25  ,45 

320 

DEG 

33  w 

LDIR  30 

340 

LABEL  -Sf'Sto* 

350 

lD  IF.  3 

350 

VIELPORT  30,115 

370 

U INCOW  -23.20,.; 

330 

frame 

393 

axes  5  , .  1  ,-zo . 

403 

AXES  5  , .  1  ,20 , 1 

410 

CSIZE  2.0. .65 

420 

MOVE  -17. 5,. 55 

430 

LABEL  ">•  <  M  }“ 

440 

A- .  62 

450 

FOR  ? - 1  TO  S 

463 

D»A-K * . 34 

470 

C$(  1  .>=  ■»:<  =  ! .  IS" 

4  SO 

CS ;  2  ■*  - -*>  1  .25" 

433 

C3( 3  >*“ *X* i . 4C" 

500 

C  5  (  4  *  “  *  X  «  !  .60" 

510 

C*( 5  .80" 

523 

CS(  6  /-"*'/  =  2 .33" 

530 

MOVE  -13.  .0 

540 

LAE EL  USING  “D  5 

53  e 

PENUP 

NEXT  K 

573 

CSIZE  1 . 

3,. 72 

s=e 

FOR  J=*l 

TO  5 

EfiS 

M* J »2 1 -2 

0 

□  02 

FDR  I=M 

TO  J *2 1 

El  2 

6< I  )-J 

520 

mode  z; I 

S  t  r  •;  I 

W7AWAV.'mVTAV.V.r^V.VV.V. 


APPENDIX  D. 


DATA  FILES  CATALOG 


i  \ A  DA  i  A  DlS\.  1*0  1 


0  17 

TDAFC2 

72387.1318  Y 

N 

10  m/s.  M  '  0.68  Tt  1 

018 

1DAFC2 

Y 

N 

input  data 

019 

FILDT2 

Y 

N 

film  cooling  data 

020 

STRFC2 

Y 

N 

St/Sto 

021 

STAVFC2 

Y 

N 

i.St/Sto  Uv 

\o 

rile  Name 

Dale 

1  il,cC. 

v  i  if  Lc  x 

Description  J 

1 

£n 

np i 

V  V..  L 

''ll  ? 

003 

004 

ID  ATI 

T  r\  .  -r*  t 

ST  A  VI 

HD  ATI 

71787.1702 

\] 

*  < 

NT 

N 

N 

X 

N 

N 

N 

15  m/s,  baseline  ; 

i 

baseline  Slav  ! 

baseline  local  h  J 

?. 

003 

TDAT2 

72087  1315 

N 

N 

10  m/s  base  line  \ 

?! 

006 

IDAT2 

Nr 

N 

■ 

3 

007 

STAY2 

X 

X 

baseline  Stav  ! 

j 

■)0  8 

IIDAT2 

\; 

V 

baseline  local  a  \ 

V, 

009 

STRDT2 

72187.1524 

N 

Y 

10  m/s  St/Sto 

> 

0 1 0 

STRV2 

N 

Y 

10  m/s  lSt/Stoiav  - 

V 

T  m 

Oil 

TDAT3 

N 

Y 

T(I)  ! 

v 

012 

IDAT3 

N 

Y 

input  data  ^ 

K 

0  1  j 

TDAT-i 

72187.1748 

N 

V 

i 

15m/sT‘l) 

K 

OH 

I  DAT  4 

N 

Y 

input  data 

R 

0  1 5 

STRDT 1 

N 

Y 

St/Sto  ! 

t 

016 

STARV1 

N 

Y 

< St/Sto >av 

IP,,K7Uf,U«,Y]iXyWv itV^WW-V-V.v  v  v:,-  -  v.v'.'>'JVW.’>,v^kiT.v  ->  ■>  ■-> 


nrrjr’J-J‘:J'-J‘J-S."J':v;-'rjv'd---r.--r*-ir.r;r*'.i  wvjntTT TXT  v  v  v 


198 


022 

TDAFC1 

72387.1622 

Y 

N 

15  m/s.  T«  1 ' 

> 

\ 

* 

023 

I D  AFC 

Y 

M 
*  -i 

,  input  OaLa 

024 

FILDT1 

Y 

M 

.  In.].  Ojla 

023 

STRFCl 

Y 

N 

,  St/S  to 

V 

026 

STAVFC1 

Y 

N 

.  1  St/Sto  »av 

\ 

r- 

027 

TDAFC4 

72487.1401 

Y 

N 

1 0  m/s,  repeat 

i 

02S 

1DAFC4 

Y 

X 

1 

»* 

& 

V 

6 

i;29 

*7  -  - 

.  i  w  U  l  '1 

t 

X 

030 

031 

w  1  >Vi  u  1 

ST  A  VFC4 

X  * 

*.  r 

S' 

.  '  i 

| 

032 

TDAFC3 

72487  1621 

y 

N 

15  m/s  repeat 

033 

iij.-vr03 

1 

j 

034 

FILDT3 

t  * 

1 

X 

(< 

f 

035 

STRFC3 

Y 

N 

i 

• 

036 

STAVFC3 

Y 

N 

£ 

037 

TDAFC6 

72587.1401 

Y 

N 

10  m/s.  M  -  1.26 

\.- 

038 

IDAFC6 

Y 

N 

. 

; 

039 

FI  LOT 6 

Y 

\ 

\ 

0^0 

STRFC6 

Y 

N 

.  St/S  to 

041 

STAVFC6 

Y 

N 

,  I St/Sto lav 

k. 

r. 

042 

TDAFC5 

72587.1808 

Y 

N 

15  m/s.  M  =  0.86 

f. 

043 

IDAFC5 

•* 

V 

4 

N 

j 

044 

FILDT5 

Y 

N 

K 

045 

STRFC5 

V 

L 

N 

St/Sto 

I 

046 

STAVFC5 

Y 

N 

.  « St/Sto iav 

P* 

047 

TDAFC7 

72887.1206 

Y 

N 

15  m/s.  M  =  0.86 

i 

jl 

048 

IDAFC7 

Y 

N 

,  Repeat 

049 

FILDT7 

Y 

N 

5 

050 

STRFC7 

Y 

N 

tf 

C 

s! 

051 

STAVFC7 

Y 

N 

t 

t 

052 

TO  A  FC8 

728X7  1555 

V 

N 

1 0  m/s  M  -  0  o,s 

•t 

4 

053 

ID  AFC8 

Y 

N 

Repeal 

V 

4 

034 

FILDT8 

V 

N 

4 

4 

055 

ST  RPC'S 

'1 

\  * 

.X 

St/Sln 

K 

056 

STAVFC8 

Y 

V 

.  * 

'St/Sto  iav 

* 


] 


3 


3, 


1 


M 


3 


3 


3 


•j 


P 


.  TOWHUV  W  '.v  '.*.'  W 


y*r»  t»  ^i  ■■ « ■  v.»  7^ 


057 

STAFC8 

Y 

. 

Su  numbers. 

0  s  s 

TDAFCIO 

72887.  P-12 

\: 

V 

!  ■ )  -rj  'o  \  1  -  f )  4  S 

()>■-) 

i  l)AK  10 

i 

i 

'ftPV  ,  Vi  **  ‘  -  2 

060 

FILDT10 

V 

X 

061 

STRFCV 1 0 

V 

Y 

St/Sto 

062 

STFCV 1 0 

Y 

Y 

<$L/Sl0  4v 

■  j  6  3 

i  */  •  ii  i  i 

—  n  Ct  a  “*  1  ^ '  *  c 

•  • 

.1  . .  . .  i 

064 

1DAFC1 1 

] 

i 

\  onex  ;en.  ~ 

Go  5 

F1LDT1  1 

V 

i 

( }  (S  o 

STRFCV !  1 

Y 

I 

St/SiO 

067 

ST  AFC  VI  ! 

y 

Y 

i St/Sto  ov 

068 

TDAFC13 

72  9  87  1345 

y 

•/ 

1 5  m/s,  M  * 

069 

ID  AFC  13 

Y 

Y 

'/one:;  ;en  ■?>  3  52  cm 

070 

FILDT13 

" 

Y 

Y 

0"1 

STRFCV 1 3 

Y 

Y 

072 

STAFCV13 

Y 

Y 

(]  /  ,5 

TD.vFC  1  5 

^2987  !  4  i  5 

\ r 

1 

^  * 

\ 

i  5  ill  /  S.  W  -'  /  •! 

074 

ID  AFC  15 

Y 

Y 

V  or  lex  gen  •«.  ,\.nn  cm 

075 

FILDT15 

’ 1 

Y 

Y 

076 

STRFCV 1 5 

" 

Y 

Y 

077 

STAFCV15 

Y 

Y 

078 

STAFC6 

72587.1 45  i 

Y 

V 

.  'i 

1 0  m/s.  VI - 1  26.  Si|  No 

079 

STAFC7 

72887.1206 

Y 

N 

1  5  m/s  VI — 0  so.  So  A' 1 

080 

STAFC3 

72487  1621 

Y 

N 

1  5  m/s.  M-0  47  Sij  \o 

081 

TDAFC2 1 

80187.1354 

Y 

Y 

15  m/s,  M  =  0.86 

082 

IDAFC21 

Y 

Y 

Vortex  ip>  6.06  cm 

083 

FILDT21 

** 

Y 

Y 

084 

STRFCV 21 

Y 

Y 

Sl/Sto 

085 

STFCV21 

Y 

Y 

averaged  St/Sto 

086 

TDAFC17 

80187  1507 

\r 

\ 

P 

I 

1  5  m/s.  M-  0  86 

087 

IDAFC17 

i 

V’ 

i 

Vortex  •<»  4.79  cm 

088 

FILDT17 

Y 

Y 

089 

STRFCV  17 

\  ~ 

I 

Y 

090 

STFC17 

Y 

Y 

199 


:3r:3\v:3-v  -•■ 


-a 

a 


>a 

[2» 

I 


SSi 


'AV.V.WJI  '-“rj'jv.v'-'j'  -/  wv.vvr. 


091 

TDAFC19 

o°: 

I D  A  PC  1 9 

09  j 

FILDT19 

094 

STRFCv  i’ 

098 

STFCV 1 9 

096 

ID  AFC  12 

097 

ID  AFC  I  2 

t  ;  r,r  i 

SO  187. 1440  V 

V 


80 1  $7  1 54 1  Y 


15  m/s.  \I  =  0.8o 
Vortex  •*  3  32  c:n 


10  m/s.  M 
Vortex  a> 


TDAFC18 

80487.1452  Y 

ID  AFC  18 

Y 

FILDT18 

Y 

STRFCY  18 

Y 

STFCV  1  8 

Y 

1 0  m/s,  M=  0  98 
Vortex  <a  479  cm 

St/Sto 

averaged  Si/Sto 


m 


*  l  ■ 

tdafczq 

80487,1529 

N  ’ 

X 

V* 

i 

10  m/s.  \i=0.98 

i  18 

IDAFC20 

Y 

Y 

Vortex  '*3  52  cm 

i  !  9 

FILDT20 

Y 

V 

x 

!  20 

ST RFC 20 

y 

t r 

i 

121 

STFCV20 

Y 

Y 

122 

TDAFC22 

80487.155  3 

Y 

Y 

10  m/s.  M=0.98 

*  ■>  ** 
t 

IDAFC22 

\r 

l 

v 

1 

Vortex  <o>  6  00  cm 

'  n  < 

.  •  l  *  \  k  *■  ,  *  ±m  km 

i 

T  ’ 

i 

V 

X 

iiFcv:: 

\r 

1 

\ ? 

1 

i  77 

TDAFC24 

80987  1943 

Y 

V 

X 

10  m/s.  M=0  98 

128 

:?.vFC2h 

V 

Y 

Vortex  <?»  h.7°  cm 

129 

FILDT24 

Y 

T  r 
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No. 

FILE  NAME 

DESCRIPTION 

131 

TPR01 

10  m/s,  M  =0.98.  heated  plate.  X=  1.867  m 

132 

TPR02 

Vortex  c<*  4.79  cm  ,  X  =  2.172  m 

133 

TPR03 

,  X  =  2.48  m 

1 3 1 

TPR04 

10  m/s,  M  -  0  98  unheated  plate.  X  -  1  8a"7  m 

138 

TPR05 

Vortex  <s>  4  79  cm.  .  X  =  2  1 72  m 

130 

TPR06 

.  X  =  2.48  m 

137 

TPR07 

,  X  =  1 .48  m 

138 

TPR08 

10  m/s,  M=0.98,  heated  plate.  X- 1 .48  m,  Vt.  <»  4.79  cm 

139 

TPR09 

lOm/s,  M=  98,  unheated  plate,  X=  1.48  m,  V  c*  3  52  cm 

140 

TPRO10 

,  V  (a  6.06  cm 

!  o  1 

TPR01 1 

10  m/s.  M=1.26.  unheated  pit  .  X  =  1 .48  m.  V  <a»  h.7'1  cm 

TPR012 


10  m/s,  M=  0  98,  unhealed  plate,  No  vortex.  V  I  ha  m 


V.1PER ATURE  PROFILES  DaTA  DISL  -2 

TPR013  10  m/s.  boundary  layer  with  film  ccoiing. 

TPR014  10  m/s,  boundarv  layer  heated  piate  oniv 

*  t  *  ;’r  .--.•••••  -  - . 

;  :  w  ai>  c.  Al-  i .  —  O.  t  io#  “x .  /  /  wui.  i  lc:.ui  v  v.i 

CAL  data,  at  X--  1  /18m 

V2LPR0  velocity  profiles 

PLOTYOR  vorticity  profiles 
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